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I.  INTRODUCTION 


BACKGROUND 

This  documentation  is  the  result  of  a  study  authorized  by  Section  22  of 
the  Water  Resources  Development  Act  (P.L.  93-251)  and  Section  214  of  the  Rivers 
and  Harbors  Act  of  1965  (P.L.  89-298).  The  study  was  requested  by  the 
Pennsylvania  Department  of  Environmental  Resources  (DER)  and  the  State  of  New 
York  Department  of  Environmental  Conservation  (DEC).  The  Philadelphia  District 
office  of  the  Corps  of  Engineers  is  responsible  for  the  conduct  of  the  work  and 
the  study  is  being  directed  by  a  group  consisting  of  staff  from  the  Corps  of 
Engineers,  the  Pennsylvania  DER,  the  New  York  DEC,  the  Delaware  River  Basin 
Commission  (DRBC) ,  and  the  United  States  Geological  Survey  (USGS).  The 
Philadelphia  District  Office  contracted  with  Water  Resources  Engineers  (WRE)  an 
operating  unit  of  Camp  Dresser  &  McKee  Inc.  and  the  firm  is  performing  work 
under  contract  number  DACW61-78-C-0127. 

A  specialized  daily  flow  model  was  developed  specifically  for  the 
Delaware  River  Basin.  This  model  is  to  be  used  to  determine  the  resulting 
discharges  in  the  Basin  for  proposed  operating  policies  for  future  and  existing 
reservoirs.  It  has  in  it  the  option  of  testing  modifications  to  existing  reser¬ 
voirs  such  as  changing  the  operating  rules  or  increasing  the  available  storage. 
Proposed  reservoirs  can  also  be  simulated  with  the  Daily  Model. 

This  report  is  a  combination  of  a  User's  Manual  and  a  Model 
Documentation.  For  those  wishing  to  understand  fully  and/or  modify  the  program. 
Chapter  II  entitled  Program  Descriptions  should  be  examined  thoroughly.  Those 
wanting  only  guidelines  to  run  the  model  are  directed  to  Chapter  III,  Input 
Descriptions,  and  Chapter  IV,  Output  Description. 


MODEL  REPRESENTATION 


The  Daily  Model  consists  of  a  branching  network  of  nodes  which  are 
linked  together.  The  nodes  represent  points  of  headwater  and  incremental 
inflows,  river  confluences,  reservoirs,  and  points  where  specific  low  flow 
objectives  are  to  be  met.  The  flow  regime  throughout  the  basin  is  determined  by 
accumulating  flows,  from  node  to  node,  in  an  upstream  to  downstream  menner 
accounting  for  headwater  and  incremental  inflows,  reservoir  releases, 
withdrawals  and  consumptive  losses  from  the  system.  Lag  or  route  times  have 
been  established  for  each  pair  of  nodes,  upstream  to  downstream.  Being  a  daily 
representation,  the  maximum  time  between  two  nodes  is  24  hours.  A  more  detailed 
discussion  of  the  development  of  the  route  times  can  be  found  in  the  Delaware 
River  Daily  Flow  Model  Report  1. 

The  model  is  driven  by  a  set  of  daily  incremental  inflows  at  each  node. 
It  incorporates  the  simulation  of  reservoir  water  supply  diversion  schemes  and 
release  rules  for  basic  conservation  at  the  reservoir  site  and  for  flow  augmen¬ 
tation  to  meet  critical  low  flow  objectives  at  Montague  and  Trenton,  New  Jersey. 
These  are  variable  and  can  be  changed  in  the  input  deck. 

In  the  following  paragraphs,  a  description  of  the  various  elements  of 
the  simulation  program  is  given. 


Flow  Objectives 

The  Daily  Model  can  currently  simulate  flow  objectives  at  any  two  loca¬ 
tions.  These  objectives  govern  the  release  schedule  of  the  corresponding  reser¬ 
voirs.  Presently,  the  two  objectives  in  operation  are  Montague  and  Trenton,  New 


^Development  of  a  Daily  Flow  Model  of  the  Delaware 
River  which  Incorporates  Reservoir  Systems  Analysis, 
Volume  I,  Camp  Dresser  &  McKee,  March  1981 
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Jersey.  Montague's  flow  governs  the  release  from  the  three  New  York  City  reser¬ 
voirs,  Pepacton,  Cannonsville  and  Neversink.  Trenton's  flow  governs  the 
releases  from  any  other  additional  reservoir(s) .  The  Daily  Model  is  capable  of 
operating  up  to  seven  more  reservoirs  for  Trenton's  objective.  If  either  objec¬ 
tive,  Montague  or  Trenton,  is  not  met  by  natural  inflows  and  conservation 
releases  from  the  reservoirs,  additional  releases  from  the  appropriate  reser¬ 
voirs  are  made  to  meet  the  objectives. 


The  New  York  City  reservoirs'  operations  are  based  on  the  daily  flow  at 
Montague.  If  the  target  is  not  met,  the  model  "backs-up"  three  days,  re¬ 
initializes  the  storage  and  the  flows,  and  makes  a  release  from  the  fullest 
reservoir.  Because  of  attenuation,  not  all  of  the  release  is  seen  at  Montague. 
Therefore,  an  input  option  that  specifies  a  multiplication  factor  to  increase 
the  release  is  available. 

The  reservoirs  used  to  meet  the  Trenton  objective  operate  in  a  much 
simpler  fashion.  Yesterday's  Trenton  flow  operates  today's  reservoir  releases. 
This  eliminates  any  overlapping  of  "back-ups"  that  could  occur  with  two  backup 
schemes.  The  attenuation  problem  is  also  handled  differently.  Because  of  the 
next  day  release  and  a  possible  maximum  two-day  lag,  the  release  made  three  days 
ago  is  subtracted  from  today's  Trenton  flow  to  better  estimate  the  natural 
Trenton  flow. 

Throughout  the  years  of  simulation,  the  model  can  simulate  dry  weather 
periods.  During  these  dry  periods,  flow  objectives  and  diversions  are  reduced. 
Therefore,  the  model  has  been  set  up  to  acknowledge  three  different  system 
conditions:  Normal,  Drought  Warning  and  Drought.  These  conditions  are  governed 
by  the  total  storage  level  of  the  three  New  York  City  reservoirs,  Pepacton, 
Cannonsville  and  Neversink.  Figure  1-1  gives  the  current  storage  scheme.  This 
scheme  is  a  modification  of  the  Drought  Emergency  Operating  Rules  -  Model  #2  for 


1-3 


the  New  York-Del aware  System  presented  in  the  Task  Group  Report,  DRBC  Docket  No. 
0-77-201 .  The  natural  hydrology  of  the  area  is  reflected  in  Figure  1-1.  For 
example,  during  the  late  winter/early  spring  months,  there  is  rainfall  activity 
causing  the  entire  basin  to  gain  water.  During  the  summer/early  fall  months, 
there  is  much  less  rainfall  and  the  basin  loses  water.  Therefore,  a  Normal 
Condition  in  October  and  November  would  be  a  Drought  Condition  for  April  through 
July,  as  is  represented  in  Figure  1-1. 

The  storage  level  is  the  tote*  1  storage  of  all  three  New  York  City 
reservoirs.  The  upper  part  of  the  graph  is  the  Normal  regime.  The  hatched 
middle  part  is  the  Drought  Warning  area.  The  lower  part  of  the  graph  is  the 
Drought  Condition.  Table  1-1  shows  the  effect  of  each  of  the  three  conditions  on 
the  Montague  objective  and  the  allowable  diversions  from  the  reservoirs  to  New 
York  City  for  use  as  water  supply.  While  currently  only  used  for  the  Montague 
objective,  the  Trenton  objective  can  also  be  linked  to  these  System  Conditions. 
Both  Table  1-1  and  Figure  1-1  contain  information  that  are  input  options  and  can 
therefore  be  changed. 


TABLE  1-1 

FLOW  OBJECTIVE  AND  DIVERSION  SCHEDULE 


Normal 

Drought  Warning 

Drought  Montague 

Objective  (cfs) 

1750 

1750 

1525 

NYC  Diversion  (mgd) 

800 

600 

430 

Ifask  Group  Report  DRBC  Docket  No.  D-77-20,  Appraisal  of 
Upper  Basin  Reservoir  Systems,  Drought  Emergency  Criteria 
and  Conservation  Measures,  DRBC,  March  1979. 
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Diversions 


The  diversions  to  New  York  City  come  from  only  Pepacton,  Cannonsville 
and  Neversink.  Any  other  reservoir  is  used  only  for  the  Trenton  flow  objective 
with  no  diversion.  Three  diversions  are  specified  in  Table  1-1  one  for  each  of 
the  System  Conditions.  These  three  diversions,  like  the  Montague  objectives, 
are  options  and  can  be  changed  for  each  simulation. 

Reservoirs 

The  reservoirs  are  operated  to  meet  the  various  flow  objectives  and 
diversion  schedules.  Currently,  only  the  three  New  York  City  reservoirs  can  be 
used  to  meet  the  Montague  objective.  Any  additional  reservoirs  are  used  for  the 
Trenton  objective.  These  reservoirs  can  be  incorporated  into  the  model  by 
merely  being  specified  in  the  input  data. 

All  reservoirs  are  located  at  headwaters.  They  are  each  set  up  with 
three  major  components:  a  continuous  storage,  a  special  release  node  and  an 
incremental  inflow  node.  When  the  reservoir  is  not  used,  the  incremental  inflow 
is  equal  to  the  flow  at  the  node.  When  in  operation,  the  reservoir's  normal 
spills  and  releases  will  also  be  lo  ted  at  the  incremental  inflow  node.  If  for 
some  reason  the  special  releases  made  to  augment  the  flows  at  the  objective 
nodes  should  not  be  released  with  the  normal  releases  conservation,  the  special 
release  node  is  specified  as  something  different  from  the  incremental  inflow 
node.  Such  is  the  case  at  Neversink  Reservoir.  Being  much  closer  to  Montague 
than  Pepacton  and  Cannonsville,  its  releases  to  augment  Montague's  flow  arrive  a 
full  day  earlier  than  releases  from  Pepacton  and  Cannonsville.  Therefore,  a 
special  release  node  has  been  specified  that  will  hold  Neversink's  releases  for 
one  day  before  allowing  them  to  be  routed  naturally  down  to  Montague.  While 
this  option  can  be  used  for  any  reservoir,  it  is  not  advised  to  do  so.  The 
Special  release  node  is  then  usually  equal  to  the  inflow  node.  The  storage  node 
is  always  specified  as  1000  plus  the  node  number  linked  to  the  incremental 
inflow  to  the  reservoir.  This  is  essential  to  keep  the  natural  inflows  and  the 
new  reservoir  releases  occurring  at  the  same  location. 
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ASSUMPTIONS  AND  LIMITATIONS 

The  Daily  Model  has  been  operated  on  only  two  systems,  a  CDC  CYBER  176 
at  United  Computing  Services  (UCS)  in  Dallas,  Texas,  and  a  CDC  CYBER  750  at  Sun 
Information  Services  in  Dallas,  Texas.  The  entire  50-year  simulation  requires 
70  CPU  seconds  at  UCS  and  160  CPU  seconds  at  Sun  Oil.  This  will  vary  slightly 
reservoir  can  be  designated  in  the  input  deck  without  changing  the  program. 
While  it  is  not  difficult  to  change  these,  it  does  require  knowledge  of  the 
program  and  FORTRAN. 

Chapter  II  contains  a  description  of  the  pertinent  program  features, 
theoretical  limitations  and  assumptions.  A  detailed  description  of  the  data 
requirements  is  provided  along  with  necessary  format  descriptions  in  Chapter 
III.  An  example  input  deck  is  given  and  described  at  the  end  of  the  Chapter. 
Chapter  IV  provides  a  discussion  of  the  output  for  the  example  input  deck. 


1-7 


-  ft. 


II.  PROGRAM  DESCRIPTION 


MODEL  STRUCTURE 

This  chapter  gives  an  overview  of  the  Daily  Model  program.  Figure  II-l 
represents  the  main  program  and  the  order  of  calling  the  nine  subroutines.  The 
main  program  and  each  subroutine  is  described  in  detail.  Included  in  the 
description  is  a  flow  chart  and  the  computer  listing. 

The  controlling  parameters  (number  of  years  in  simulation,  monthly 
counter,  inflow  tape  number,  etc.)  are  initialized  and  read-in  (INDATA),  the 
nodes  are  ordered  from  upstream  to  downstream  and  then  linked  to  their  incremen¬ 
tal  inflow  (ORDER  and  LINKER),  all  before  RGFL  (main  program)  enters  the  daily 
and  yearly  loops.  The  System  Condition  storage  levels  are  also  calculated 
( INTRPL)  before  entering  the  loop.  The  incremental  inflows  are  read  for  the 
current  day  (INFLOW).  The  upstream  flows  are  added  to  the  downstream  incremen¬ 
tal  inflows  (ROUTE)  once  the  reservoirs  have  been  operated  (RESOP).  If  Montague 
needs  more  water,  the  simulation  backs  up  three  days  (NEWFLO).  After  each  day 
which  is  not  in  the  back  up  procedure,  all  the  flows  are  written  on  a  special 
tape,  later  to  be  reduced  into  tabular  form  (OUTPUT  or  MONTH). 


RGFL  Program 

The  main  program  of  the  Daily  Flow  Model  is  called  RGFL,  an  acronym  for 
the  Regulated  Flow  model.  It  calls  most  of  the  subroutines  and  keeps  track  of 
the  daily  and  yearly  counters,  as  well  as  the  backup  procedure.  Figure  1 1 -1 
gives  a  general  flow  chart  for  the  RGFL  program.  Parameters  are  initialized  and 
variables  are  read  in  (INDATA)  before  the  model  begins  the  daily  loop.  The 
subroutine  LINKER  connects  the  model  nodes  to  the  incremental  inflow  nodes  on 
the  tape. 


II-l 


Once  inside  the  daily  loop,  RGFL  calls  INFLOW  which  reads  the  incremen¬ 
tal  inflows  and  ROUTE  which  moves  the  water  through  the  basin  with  lag  functions 
and  various  reservoir  operations.  The  past  three  days'  flows,  storages  and 
counters  are  then  set  to  be  used  if  and  when  backup  occurs. 

The  main  purpose  of  RGFL  is  to  check  whether  or  not  the  model  needs  to 
make  a  release  to  augment  the  Montague  flow.  If  a  release  is  required,  RGFL 
calls  NEWFLO  which  backs  up  the  simulation  three  days.  A  flag  is  set  in  NEWFLO 
such  that  the  model  can  never  backup  to  that  particular  day  again.  This  flag  is 
cnecked  in  RGFL  after  the  past  three  days'  flow  values  have  been  set.  In  the 
flowchart,  this  check  occurs  in  two  places.  The  question  "Is  the  Model  in 
Backup?"  directly  checks  the  flag.  The  next  one  "Should  the  Model  Backup?"  in 
actuality  is  two  questions.  First,  it  asks  if  the  Montague  flow  is  lower  than 
the  objective.  If  it  is,  it  asks  if  the  model  just  finished  backing  up.  In 
other  words,  whether  the  simulation  is  at  the  point  where  it  was  required  to 
backup  before.  If  so,  the  model  will  not  backup  due  to  the  check  on  the  backup 
flag  and  instead  writes  on  a  scratch  tape  all  the  flows,  storages  and  special 
outputs  from  three  days  ago.  By  writing  the  information  from  three  days  ago, 
there  is  no  possible  way  to  rewrite  days  due  to  the  backup  procedure. 


If  the  model  beginning  the  backup  procedure,  it  goes  back  to  INFLOW  in 
order  to  set  the  inflows  for  the  current  day.  This  is  done  without  incrementing 
the  daily,  monthly  or  yearly  counters.  It  then  calls  ROUTE  which  makes  the 
release  from  the  fullest  NYC  reservoir  on  the  first  day  only. 

RGFL  performs  two  additional  calculations.  The  first  calculation  is 
the  release  needed  to  maintain  the  Trenton  flow  objective.  The  second  calcula¬ 
tion,  when  appropriate,  is  Merrill  Creek's  operation.  Designed  in  reality  to 
release  whenever  the  Trenton  flow  is  less  than  3100  cfs,  RGFL  will  increase 
Merrill  Creek's  target  of  3100  cfs  if  the  Trenton  objective  is  greater  than  3100 
cfs.  This  is  discussed  in  much  greater  detail  in  the  Delaware  River  Daily  Flow 
Model  Report  1. 


^Development  of  a  Daily  Flow  Model  of  the  Delaware 
Riyerwhich  Incorporates  Reservoir  Systems  Analysis, 
Volume  I,  Camp  Dresser  &  McKee,  March  1 98l 


The  daily  and  yearly  loops  are  continued  until  the  number  of  years  in 
the  simulation  have  been  completed  or  the  end  of  the  inflow  tape  has  been 
reached.  At  this  point,  either  MONTH  or  OUTPUT  is  called  to  create  the  tabular 
output.  MONTH  produces  the  monthly  average  values  while  OUTPUT  produces  the 
daily  values  for  each  year. 
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Figure  I I - 1 .  Flowchart  for  Main  Program  RGFL. 
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REGULATED  FLOW  MODEL 

DAILY  FLOW  MODEL  OF  THE  DELAWARE  RIVER  BASIN 
WRITTEN  BY  J.W.RIDGUAY 

REVISED  BY  K.E. ZIMMERMAN  DECEMBER  1979 

REVISED  BY  S. A. HANSON  FEBRUARY  1981 

********************************************************************** 

PROGRAM  RGFL( INPUT » OUTPUT »TAPE5= INPUT » TAPE6=0UTPUT»TAPE10*TAPE20i 
.  TAPE30»TAPE99  > 

********************************************************************** 

TAPE  CLARIFICATION J 

TAPE10*  BINARY  OUTPUT  OF  SIMULATION 
TAPE20I  EXTRA  TAPE  -  NOT  CURRENTLY  USED 
TAPE305  INCREMENTAL  INFLOW  TAPE  -  REQUIRED 
TAPE99J  EXTRA  TAPE  -  NOT  CURRENTLY  USED 
********************************************************************** 
COMMON  /CONTROL/  NYEAR, NTAPE rNOUT » NOG .NNODE. ITAPE f I YRf IMO. IDY 
COMMON  /CHANL/  NGT0C100) fNAHEG< 100) .NGTOG( 100.10) fNOUTPTOOO) f 

.  LINK ( 100) fG( 100) fQIN( 100) »QFK( 4) fQFK0(4) fQPR(4) f 

,  OPRO ( 4 ) 

COMMON  /RESERV/  NRESfXFLOWf IFLAGfVOL( JO) fLAGRES(IO) fLNKRES(IO) f 

,  EVAP< 12f 10 ) f STORE ( 10 ) >  CONS (12f10)fQ01(100)f 

.  002(100) fNOGG f  ICK  f  FLOWL  fZFLOWf ICK1 

COMMON  /STOR/  STOREl ( 10) fST0RE2< 10) »ST0RE3< 10) fST0RE4( 10) 

COMMON  /LAG/  OOUT ( 100) fXM< 1 00) fB( 100) 

COMMON  /DATE/  I YR3* IM03» IDY3f IYR2f IM02r IDY2f IYR1 f IM01 f IDY1 f 
.  01 ( 100) fQ2( 100) 

COMMON  /PHAS2/  NYCRESf ICHK»ST0REN(31 » 12) f STORED (31 f  12)  f 

.  IHVERS<3>f0BJECT(3)fEVF(12f10)fIEVPfDFL0W(3>f 

.  TARG<3) 

COMMON  /PADER/  ACONS< 12* 10 )» TRENT » YFLOWf REL ( 1 0) fNSTEPfMDATE 
COMMON  /FLAGS/  NTRfNDRfMERfMPRfNTGfMRFSfNBELfMBL 

COMMON  /CAL/  CALIB1 , JOUTf JMIN 

COMMON  /OUTS/  NSPO 

COMMON  /MOMIN/  QMM < 1 2 f50f 10 > fQMIN( 10) 

COMMON  /BACKUP/  YF5»YF4fYF3fYF2fYF1 

DATA  QM IN/ 101999999.  / 

C  *****  INITIALIZE  PARAMETERS 
IFLAG=2 
IEVP=1 
ICHK=1 
FL0WL=0 . 

XFL0W=0. 

YFL0W=0. 

ZFLOW=0. 

NSTEP^O 

MDATE=0 

C  *****  CALL  INDATA  (READS  INPUT  DATA) 

CALL  INDATA 
IM0=1 
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idy=31 

c  *****  CALL  LINKER  (ASSIGNS  INCREMENTAL  INFLOWS  TO  NODES) 

CALL  LINKER 
NOGG=NOG 

NOG=NOG+NRES+NSPO 

C  *****  NT APE  IS  AN  UNFORMATTED  OUTPUT  TAPE 
WRITE(NTAPE)NOG 
WRITE(NTAPE) (NAMEG( I ) » 1=1 .NOG) 

C  *****  ENTER  DAILY  LOOP  FOR  NYEAR  YEARS 
DO  1000  1=1. NYEAR 
DO  800  11=1.366 

C  *****  CALL  INFLOW  (READS  INFLOWS  AT  NODES  AND  RESERVOIRS.  CALCULATES 

C  RESERVOIR  LEVELS) 

10  CALL  INFLOW 

C  *****  CALL  ROUTE  (SETS  TARGETS.  OPERATES  RESERVOIRS  AND  ROUTES  FLOWS 

C  DOWNSTREAM) 

CALL  ROUTE 

C  *****  SAVE  THE  PAST  THREE  DAYS  DATES 
IYR3=IVR2 
I YR2=I YR1 
I YR1=IYR 
IMD3=IM02 
IM02=IM01 
1MO1=IM0 
IDY3=IDY2 
IIiY2=IDYl 
IDYJ  =IDY 

lF(NRES.EQ.O)  GO  TO  610 

C  *****  SAVE  THE  PAST  THREE  DAYS  STORAGE  LEVELS 
DO  600  N=1 . NRES 
ST0RE4(N)=ST0RE3(N) 

ST0RE3(N)=ST0RE2(N) 

ST0RE2(N)=ST0RE1(N) 

ST0RE1 (N)=STORE(N) 

600  CONTINUE 
610  CONTINUE 

C  *****  SAVE  THE  FAST  THREE  DAYS  FLOWS 
DO  700  N=1 .NOG 
002( N ) =001 (N) 

001 (N)=02(N> 

Q2(N)=Q1(N) 

Q1(N)=Q(N) 

700  CONTINUE 

C  *********************************************************************** 

C  LOWER  BASIN  TARGET  IS  AT  NODE  ICK1 

C  TRENT  IS  THE  TARGET  AT  2ND  TARGET  NODE  (ORIG.  TRENTON) 

C  YFLOW  IS  THE  AMOUNT  THE  2ND  TARGET  NODE  MISSED  ITS  FLOW 

C  ZFLOW  IS  FOR  THE  MERRILL  CREFK  OPTION 

C  *********************************************************************** 

YF6=YF5 

YF5=YF4 

YF4=YF3 

YF3«YF2 

YF2=YF1 
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YF1=YFL0W 

YFLOW= < TRENT-01 (NTR ) )  +YF3 
IF(YFLOW)  701»701»702 

701  YFL0W=0. 

702  CONTINUE 

C  *****  CHECK  IF  TARGET  FLOW  IS  ACHIEVED  AND  IF  BACKUP  HASNT  OCCURRED 
C  *****  NOTE  :  THIS  OCCURS  ONLY  WHEN  NYC  RESERVOIRS  ARE  INCLUDED 
IF<  Q< ICK) .LT . FLOWL . AND . IFLAG.LT .0)CALL  NEWFLO 
IF< IFLAG.GT .0)  GO  TO  695 

WRITE (NT APE)  IYR3»IH03»IDY3»<Q01(N> .N=1»N0G> 

695  CONTINUE 

C  *****  IFLAG  KEEPS  TRACK  OF  THE  NUHBER  OF  DAYS  SINCE  THE 
C  *****  LAST  BACKUP  OCCURRED 
IFLAG=IFLAG-1 
ZFL0W=Q1 (NTR)-3100 

C  *****  NEXT  STATEHENT  ALLOWS  FOR  A  VARIABLE  HERRILL  CREEK  OPTION 
IF( TRENT *GT >3100)  ZFLOW  =Q1 (NTR)-TRENT 
IF ( I FLAG) 800 r 10 » 10 
800  CONTINUE 
1000  CONTINUE 

IF (JOUT >1100(1200 >1300 
1100  CALL  OUTPUT 
GO  TO  1400 
1200  CALL  MONTH 
GO  TO  1400 
1300  CALL  MONTH 
CALL  OUTPUT 
1400  CONTINUE 
STOP 
END 
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Subroutine  INDATA 


The  main  function  of  INDATA  is  to  read  and  write  all  of  the  card  input 
variables  as  seen  in  Figure  1 1 -2 .  If  INDATA  detects  any  inconsistencies  in  card 
order,  an  error  message  is  printed  and  execution  is  immediately  terminated.  It 
is  therefore  very  important  to  put  in  the  flags  specified  in  the  various  data 
series.  INDATA  calls  the  subroutine  ORDER  which  connects  the  upstream  to 
downstream  nodes.  The  reservoir  evaporation  is  changed  from  inches/month  to 
cfs-days  (a  volume)  right  after  the  reservoir  information  is  read  in.  INDATA 
also  calls  INTRPL  which  calculates  the  reservoir  levels  for  the  System 
Conditions.  INDATA  returns  to  RGFL. 
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Figure  1 1 -2 -  Flowchart  for  Subroutine  INDATA 
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SUBROUTINE  INDATA 

c  m*******************************?************************************ 

C  INDATA  READS  THE  INPUT  DATA' 

C  YEARS  OF  SIMULATION 

C  NUMBER  OF  RESERVIORS 

C  MODEL  NODES  AND  LAG  FUNCTION  COEFFICIENTS 

C  OUTPUT  OPTIONS 

COMMON  /CONTROL/  NYEARfNTAPE  fNOUT fNOGfNNODEf I  TAPE  f I YRf IMOf IDY 
COMMON  /CHANL/  NGTO ( 1 00 ) . NAMEG( 100) fNGTOG< 100 f 10) f NOUTPT (100 ) , 

.  LINK ( 100 ) fQ( 100 ) f QIN( 1 00) fQFK ( 4 ) f  QFKO ( 4 ) f  QPR ( 4 ) f 

.  QPR0<4) 

COMMON  /RESERV/  NRESf XFLOWf IFLAGf VOL ( 1 0 ) f LAGRES! 10 ) f LNKRES ( 10 ) f 
.  EVAP( 12f 10) f STORE! 10)  f CONS! 12»10) rQ01< 100) » 

.  Q02( 100) iNOGGfICKfFLOWLfZFLQWf ICK1 

COMMON  /STOR/  ST0RE1 < 10 > f ST0RE2! 10 > f ST0RE3! 10) f ST0RE4 ( 10 ) 

COMMON  /INOUT/  NTITLE1 ( 10) fNTITLE2( 10) fNTITLE3( 10 f 100) 

COMMON  /LAG/  QOUT < 100 ) f XH! 100) f B ( 100 > 

COMMON  /CAL/  CALIB1 f JOUTf JMIN 

COMMON  /PHAS2/  NYCRESf ICHK fSTOREN! 31 f 12) fSTORED! 31 f 1 2 ) f 
.  DIVERS!  3 ) f OBJECT (3>fEVP(12f10)f IEVPf  DFLOW ( 3 ) » 

.  TARG<3) 

COMMON  /PADER/  ACONS< 12 . 10) , TRENT f YFLOWfREL ( 10) f NSTEPfMDATE 
COMMON  /CONSUM/  NCONf  ICON!  i»0>  f  CLQSS!  10f  12  > 

COMMON  /OUTS/  NSPO 
COMMON  /WALLIN/  MWL 

COMMON  /FLAGS/  NTRfNDRfMERfMPRfMTGfMRESfNFEL  fMBL 

DIMENSION  MM0N(12)fNAMR(5) 

DATA  MM0N/31f28f31f30f31f30f31f31f30f31f30.31/ 

WRITE(6f 10) 

10  FORMAT ( 1H1 ) 

C  *****  READ  TWO  TITLE  CARDS  A.K.A.  ’A'  CARDS 

READ( 5  f 1300) AC  fNTITLEI 
IF(AC.NE.IHA)  GO  TO  5100 
WRITE( 6 » 1301 ) AC  fNTI TLE 1 

1300  FORMAT (A1f4Xf10A4) 

1301  FORMAT (1XfA1f4Xf10A4) 

READ(5f 1300)  ACfNTITLE2 
IF(AC.NE.IHA)  GO  TO  5100 
WRITE(6f1301)ACfNTITLE2 
ICK=0 

ICK1=0 

C  *****  READ  ONE  CONTROL  CARD  A.K.A.  'B'  CARD 

READ(5f1000)BCf ITAPEfNTAPEfNYEARfCALIBI fNYCRESf JOUT  f JMIN f MWL 
IF(BC.NE.IHB)  GO  TO  5200 
IFCCALIBl.EQ.O.)  CALIB1=1.25 
IF(NYEAR.EQ.O)  NYEAR=1 

WRITE(6f299B)BCfITAPEfNTAPEfNYEARfCALIB1fNYCRESfJ0UTfJMINfMWL 
1000  FORMAT! A1 f4Xf3I5fF5,2f4I5) 

2998  F0RMAT(1XfA1f4Xf3I5fF5.2f4I5) 

C  *****  READ  UP  TO  100  CHANNEL  CARDS  A.K.A.  'C*  CARDS 

DO  100  1=1 f 100 
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READ(5r 1100)CCrNAMEG( I ) *NGTO( I )  r XM( I ) rB< I ) r ICKK 
IF(CC.NE.IHC)  GO  TO  150 

WRITE (&r 1101 )CCrNAMEG< I ) ?NGTO( I )>XN(I)iB(I) » ICKK 
IF(NAMEGd) )5800»1120» 1110 
1100  F0RMAT<Alr4Xr2I10rE10.4rF10. Or  110) 

1101  F0RMATdXrAlr4Xr2I10rE10.4rF10.0rI10) 

1110  IFdCKK.EQ.2)  ICK1=I 
IF( ICKK ,EQ. 1 ) ICK=I 
Bd  )=B( I )/24 . 

XMd  )  =  XMd )/24 . 

100  CONTINUE 
150  CONTINUE 
1120  N0G=I-1 

CALL  ORDER 
IFdCK.EO.O)ICK=NOG 
IFdCKl.EG.O) ICK1=N0G 
DO  2000  N=  1  r 100 

C  *****  READ  UP  TO  100  CHANNEL  OUTPUT  CARDS  A.K.A.  *D*  CARDS 

C  *****  NODES  FOR  WHICH  YOU  REQUEST  OUTPUT 

READ(5rl200)DCrN0UTPT(N)r(NTITLE3(NNrN)rNN=lrl0) 

1200  FORMAT (Alr4XrI5rl0A4) 

IF < DC . NE . 1HD )  GO  TO  5810 

WRITE(Arl201)DCrN0UTPT<N)r(NTITLE3(NNrN)rNN=lrl0) 

IF (NOUTPT (N) )5810r2100r 2000 
1201  FORMAT (lXrAlr4XrI5r 10A4 ) 

2000  CONTINUE 
2100  CONTINUE 
NOUT =N-1 

C  *****  ENTER  INFORMATION  FOR  UP  TO  10  RESERVOIRS  A.K.A  ’E'  CARDS 
DO  3000  N=1 » 1 1 
NX=N+NOG 

READ(5r2200)ElCr IRESr VOLUME rLRESr  AREA r (NAMR( IR) * 1R=1 r5  > 

IF ( E1C . NE . 2HE 1 )  GO  TO  5820 

WRITE ( 6r2201)ElCr IRESr VOLUME rLRESr AREAr ( NAMR< IR > r IR=1 r 5 ) 

IF ( IRES)  5820r3500r2220 
2200  FORMAT (A2r3XrI5rE10.4rI10rF10,0r5A4) 

2201  FORMAT (lXrA2r3XrI5rE10,4rI10rF10.0r5A4) 

2220  IF(IRES.EQ.2155)  MRES=N 

C  *****  FIND  THE  NODE  WHERE  THE  REQUIRED  RESERVOIR  RELEASES  ARE  MADE 
DO  2210  NN=1 rNOG 
IF (NAMEG(NN ).EQ.LRES)GO  TO  2215 

2210  CONTINUE 

WRITE (6r 221 1)LRES 

2211  FORMAT ( 19H  NO  MATCH  FOR  LRES  rI5> 

STOP 

2215  LAGRES(N)=NN 

IFILRES.EQ. 1050>MPR=NN 
IF < LRES. EQ. 1 4501 MBL=NN 
VOL ( N ) =V0LUME/86400 . 

STORE<N)=VOL(N) 

STORE 1 (N)=STORE(N> 

ST0RE2(N)=ST0RE(N) 

ST0RE3(N)=ST0RE(N> 
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ST0R£4(N)=STGRE(N) 

C  *****  READ  3  OTHER  RESERVOIR  CARDS 
C  *****  FIRST  IS  THE  BASIC  CONSERVATION  RELEASE 
READ(5i2250) <CONS( I»N)*I=1»12) 

WRIT£(6 . 2255  > (  CONS< I»N)»I=1»12) 

C  *****  SECOND  IS  THE  AUGMENTED  CONSERVATION  RELEASE 
READ<5»2250) (ACONS( I.N). 1=1.12) 

WRITE(6»2255) ( ACONS( I.N). 1=1.12) 

C  *****  THIRD  IS  THE  EVAPORATION  IN  INCHES  PER  MONTH 
READ (5. 2250) (EVAP( I.N) .1=1.12) 

WRITE < 6.2255 >(EVAP( I.N) .1=1.12) 

C  *****  CHANGE  EVAPORATION  FROM  INCHES/MONTH  TO  CFS-DAYS 
C  *****  0.042=43560/(12*86400) 

DO  2260  1=1.12 
AMON=MMON( I > 

2260  EVAP(  I » N) =EVAP( I »N)*AREA*0 . 042/AM0N 
WRITE <6.2 255)(EVAP(I»N).I=1*12) 

C  *****  STORE  THE  EVAPORATION  VALUES  TO  BE  USED  WITH  NYC  INFLOWS 
DO  2245  IM= 1 . 12 
EVP(IM,N)=EVAP(IM.N) 

2245  CONTINUE 

2250  FORMAT < 5X » 12F5 . 0 ) 

2255  FORMAT < 12F10. 2 ) 

IRES=IRES-1000 

C  *****  FIND  THE  NODE  WHERE  CONSERVATION  RELEASES  AND  OVERFLOWS  ARE  MADE 
DO  2300  NN=1 .NOG 
IF (NAMEG(NN) .EG. I RES) GO  TO  2350 
2300  CONTINUE 

IRES=IRES+1000 

WRITE(6.2330)IRES 

2330  FORMAT ( 23H  NO  MATCH  FOR  RESERVOIR  .15) 

STOP 

2350  LNKRES(N)=NN 

NAMEG<  NX )  =  IRES+ 1000 
3000  CONTINUE 

WRITE(6.3100) 

3100  FORMAT (21H  TOO  MANY  RESERVOIRS  ) 

STOP 

3500  CONTINUE 
NRES=N-1 

IF ( NYCRES.LT  .  1 )  GO  TO  3615 

C  *****  CALL  INTRPL  WHICH  CALCULATES  THE  NORMAL  AND  DROUGHT  CURVES 
C  *****  FOR  THE  NYC  RESERVOIRS 
CALL  INTRPL  (STOREN) 

CALL  INTRPL  (STORED) 

C  *****  READ  IN  DIVERSION  A.K.A.  'H*  CARD 

READ! 5. 3600)  CH. (DIVERS( I ) . 1=1 .3) 

IF(CH.NE.IHH)  GO  TO  5300 

WRITE (6. 3601 >CH .(DIVERS (I). 1=1.3) 

3600  FORMAT (A1.4X.6F10.0) 

3601  FORMAT (1X.A1.4X.6F10.0) 

C  ******  READ  IN  OBJECTIVE  CARD  A.K.A.  *1*  CARD 

READ  ( 5. 3600) Cl » ( OBJECT < I ) » 1  =  1 .3) 
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IF(CI.NE.IHI)  GO  TO  5400 
WRIT£(6, 3601 )CI , (OBJECT ( I ) *  1  =  1  *  3) 

DO  3610  K=1 »3 

C  *****  CONVERT  DIVERSIONS  FROM  HGD  TO  CFS 
DIVERS! K)=DIVERS!K> *1.5473 
3610  CONTINUE 
3615  CONTINUE 

C  *****  READ  THE  SPECIAL  OUTPUT  SELECTIONS  A.K.A.  ’J1  CARDS 

READ!5,3999>CJ,NSPQ 
3999  FORMAT (AS, 15) 

WRITE! 6, 3979) CJ,NSPO 
3979  F0RMAT(1X»A5,I5) 

IF(CJ.NE.IHJ)  GO  TO  5900 
IF  <  NSPO >4011,4011, 3989 
3989  DO  4001  N=1,NSP0 
NNG-N064NRES4N 
N0UT=N0UT+1 

READ(5* 1200)  CJ,NOUTPT (NOUT ) * (NTITLE3(NNrN0UT ) »NN=1 » 10) 

WRITE (6, 1201 )C J  »NOUTPT ( NOUT ) » ( NTITLE3! NN,NOUT ) » NN=1 , 10 ) 

NAMEG ( NNG ) =NOUTPT ( NOUT ) 

4001  CONTINUE 

C  *****  READ  DIR  CANAL  FLOW  CHANGES  A.K.A  *K‘  CARD 

4011  READ< 5, 1007 )CK » ( DFLOW! I > , 1=1,3) 

IF(CK.NE.IHK)  GO  TO  5500 
WRITE(6,1027)CK,< DFLOW! I >,1=1,3) 

1007  FORMAT  <A1»4X»3F10.0) 

1027  FORMAT (1X,A1,4X,3F10.0) 

IF(NYCRES)  5010,5010,5020 

C  *****  READ  TRENTON  TARGETS  FOR  OPERATING  RULES  A.K.A.  *L '  CARD 
5020  READ<5, 1007 ) CL , ( TARG! I) ,1=1*3) 

IF(CL.NE.IHL)  GO  TO  5600 
WRITE! 6 , 1027) CL , ( TARG ( I ) , 1  =  1 , 3 ) 

C  *********************************************************************** 

C  READ  INPUT  FOR  CONSUMPTIVE  LOSSES  A.K.A.  'M'  CARDS 

C  NCON  =  NUMBER  OF  NODES  WITH  CONSUMPTIVE  LOSSES 

C  ICON  =  NODE  NUMBER 

C  CLOSS  =  LOSSES  SPECIFIED  FOR  EACH  MONTH  FOR  EACH  NODE 

C  *********************************************************************** 
5010  READ<5» 2001 )CM, NCON 
IF (EOF! 5) )2003» 1999 
1999  CONTINUE 

IF(CM.NE.IHM)  GO  TO  5700 
WRITE <6, 9999) CM, NCON 

2001  F0RMAT(A1 »4X» 110) 

9999  FORMAT (1X,A1»4X*I10) 

DO  6000  1=1, NCON 

RE AD <5, 2002) ICON! I l , (CLOSS! I , J) , J=1 , 12) 

6000  CONTINUE 

2002  F0RMAT(5X, 110, 12F5.0) 

C  *****  FIND  THE  NODES  WHERE  THE  CONSUMPTIVE  LOSSES  ARE  MADE 
DO  6010  1=1, NOG 
DO  6010  K=1 ,NCON 

IF! ICON(K) .EQ. NAMED! I) )  ICON(K)=I 


6010  CONTINUE 

DO  6020  I=1fNCON 

WRITE (6> 2002) ICON < I ) r <CLOSS( I » J) » J=1 » 12) 

6020  CONTINUE 
RETURN 
2003  NC0N=0 
RETURN 

5000  CONTINUE 
RETURN 

5100  WRITEC6.5101) 

5101  FORMAT (54H  ERROR  t  TITLE  *A*  CARD  LEFT  OUT  5  PROGRAM  TERMINATED) 
STOP 

5200  WRITE(6»5201) 

5201  FORMAT (55H  ERROR  *.  CONTROL  CARD  *B*  MISSING  :  PROGRAM  TERMINATED) 
STOP 

5300  WRITE(6.5301) 

5301  FORMAT (59H  ERROR  !  NO  DI VERSIONS-CHECK  ,Hi  CARD  5  PROGRAM  TERMINA 
.TED) 

STOP 

5400  URITE<6*5401 ) 

5401  FORMAT (54H  ERROR  5  NO  OBJECTIVE  *1*  CARD  \  PROGRAM  TERMINATED  ) 
STOP 

5500  WRITE ( 6»5501 ) 

5501  F0RMAT(57H  ERROR  1  NO  DIVERSIONS  TO  DIR  CANAL  5  PROGRAM  TERMINATE 
.D  ) 

STOP 

5600  WRITE(6> 5601 ) 

5601  FORMAT (68H  ERROR  5  NO  2ND  TARGET  FLOWS  -  CHECK  'L*  CARD  J  PROGRAM 
.TERMINATED  ) 

STOP 

5700  WRITE(6»5701) 

5701  FORMAT (69H  ERROR  J  INCORRECT  INPUT  FOR  CONSUMPTIVE  LOSSES  5  PROGR 
.AM  TERMINATED) 

STOP 

5800  WRITE( 6  »5801 ) 

5801  FORMAT (57H  ERROR  J  NO  ’ZERO*  TO  END  'C*  CARDS  *.  PROGRAM  TERMINATE 
.D  ) 

STOP 

5810  WRITE (6  r 581 1 ) 

5811  FORMAT (63H  ERROR  i  INCORRECT  FORMAT  OF  ‘D1  CARDS  l  PROGRAM  TERMIN 
.ATED  ) 

STOP 

5820  WRITE(6»5821 ) 

5821  FORMAT (77H  ERROR  *.  INCORRECT  FORMAT  OF  RESERVOIR  INPUT  'E*  CARD  *. 

.  PROGRAM  TERMINATED  ) 

STOP 

5900  WRITE<6.5901> 

5901  FORMAT (53H  ERROR  5  SPECIAL  OUTPUT  MISSING  J  PROGRAM  TERMINATED  ) 
STOP 

END 

9 
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Subroutine  ORDER 


Subroutine  ORDER'S  flow  chart  is  shown  in  Figure  1 1-3 .  ORDER  connects 
the  upstream  to  downstream  nodes  by  use  of  the  order  of  the  nodes,  or  "slots" 
that  they  fill.  Specifically,  a  zero  array  is  filled  in  such  that  all  model 
nodes  are  connected  to  all  those  nodes  flowing  into  them. 

The  process  of  filling  in  the  array  is  geared  to  the  downstream  node. 
Each  downstream  node  is  checked  against  all  the  model  nodes.  When  a  match  is 
made,  the  2-dimensional  zero  array  sets  the  appropriate  element  to  a  positive 
integer.  The  first  subscript  is  the  matched  model  node  slot  number.  The  second 
subscript  is  the  counter  of  the  upstream  nodes.  The  element  is  equal  to  the 
slot  number  of  the  upstream  node  corresponding  to  that  particular  downstream 
node  being  checked.  This  results  in  a  2-dimensional  array  where  the  first 
subscript  corresponds  to  a  model  node,  the  second  subscript  to  the  number  of 
nodes  flowing  into  it  and  the  elements  equaling  those  particular  upstream  nodes. 
Therefore,  a  headwater  node  would  have  only  zeroes  as  its  elements. 

The  following  is  a  simplified  example  of  the  process  of  ORDER.  Given  a 
reach  with  six  nodes  as  in  the  network  shown  in  Figure  II-3A,  the  upstream  to 
downstream  data  would  be  input  as  in  Table  II-l. 


TABLE  II-l 
SAMPLE  INPUT 


Upstream  Downstream 


1 

2 

3 

4 

5 

6 


3 

3 

5 

5 

6 
0 


Pair  Number 

(Pair  1) 
(Pair  2) 
(Pair  3) 
(Pair  4) 
(Pair  5) 
(Pair  6) 


Figure  II-3A.  Sample  Network 
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When  ordering  the  nodes,  the  process  examines  each  pair  of  nodes.  In 
the  first  pair,  1  is  the  upstream  node  and  3  is  the  downstream  node.  3  is  com¬ 
pared  to  all  the  model  (upstream)  nodes.  It  is  equal  to  the  upstream  node  in 
the  3rd  pair  in  the  list.  The  first  element  of  the  array  to  be  filled  in  will 

be:  (3,1)=1,  which  says  that  the  first  node  to  flow  into  3  is  1.  The  second 

pair  of  nodes  is  then  examined.  Again  the  downstream  node  is  3  while  the 
upstream  node  is  now  2.  3  is  again  compared  to  all  the  model  nodes.  Again  it 
matches  the  third  pair.  The  next  element  to  be  filled  in  will  be  (3,2)=2.  This 

says  that  the  second  node  to  flow  into  3  is  2.  The  third  pair  is  examined.  The 

upstream  node  is  now  3  and  the  downstream  node  is  5.  5  is  compared  to  the 

model  nodes  and  matches  pair  number  5.  The  third  element  to  be  filled  in  will 

be  (5,1 )=3.  The  first  node  to  flow  into  5  is  3.  The  fourth  pair  is  examined,  5 

is  compared  and  matched  to  the  fifth  pair.  The  fourth  element  to  be  filled  in 

is  (5,2)=4.  The  second  node  to  flow  into  5  is  4.  The  last  full  pair  is  exa¬ 

mined.  The  last  element  to  be  filled  in  is  (6,1 )=5.  All  elements  not  being  set 
specifically  to  a  positive  value  are  equal  to  zero.  The  final  array  is  given  in 
Table  1 1-2.  ORDER  returns  to  INDATA. 

TABLE  1 1 -2 
SAMPLE  MATRIX 

Position  of  Nodes  Flowing  into  Model  Node 
1  2  3  4  5 

o  0  0  0  0 

0  0  0  0  0 

1  2  0  0  0 

0  0  0  0  0 

3  4  0  0  0 

5  0  0  0  0 


Model  Node 

- T - 

2 

3 

4 

5 

6 
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Figure  1 1-3.  Flowchart  for  Subroutine  ORDER. 
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SUBROUTINE  ORDER 

c  *********************************************************************** 

C  ORDER  CREATES  AN  ARRAY  THAT  LINKS  DOWNSTREAM  NODES  TO  THE 

C  UPSTREAM  NODES  FLOWING  INTO  IT. 

c  *********************************************************************** 

C  *****  ORDER  CREATES  AN  ARRAY  WITH  THE  NODES  UPSTREAM  OF  EACH  NODE 
COMMON  /CONTROL/  NYEAR,NTAPE.NOUT,NOG,NNQDE»ITAPE»IYR#IMO>IDY 
COMMON  /CHANL/  NGTO< 100) .NAMEG! 100) .NGTOG< 100. 10) .NOUTPT! 100) . 

.  LINK(100).0(100)>OIN(100)»QFK(4) »QFKO( 4) »QPR( 4) * 

.  QPRO ( 4 ) 

NG=100 

NLETS=0 

DO  750  N=1 t NOG 
NN=NOG+NLETS 

C  *****  FIND  WHICH  NODE  EACH  NODE  IS  UPSTREAM  OF  AND 
C  *****  STORE  IN  THE  NGTOG  ARRAY 
DO  720  NG0T0=1 *  NN 

IF  (NGTO(N) .EQ.NAMEG(NGOTO) )  GO  TO  730 
720  CONTINUE 

C  *****  CREATE  DUMMY  GUTTERS  AS  NEEDED 

C  *****  NLETS  IS  THE  NUMBER  OF  END  NODES  IN  THE  MODEL  RUN 
NLETS=NL£TSM 
NGOTO=NOG+NL£TS 
NAMEG !NGOTO)=NGTO(N) 

IF  (NGOTO.GT.NG)  GO  TO  760 
730  CONTINUE 

DO  740  J=1 . 10 

C  *****  BUILD  THE  NGTOG  ARRAY 

IF ( NGTOG (NGOTOrJ) . GT.O)  GO  TO  740 
NGTOG (NGOTO«J)=N 
GO  TO  750 
740  CONTINUE 
750  CONTINUE 

NOC=NOG+NLETS 
GO  TO  780 

C  *****  ERROR  IN  DATA 

7A0  URITE<6»770)  NGOTOr NG 

770  FORMAT (45H  ****  ERROR  ****  THE  ASSIGNED  CHANNEL  NUMBERS  »I5» 

1  24H  WHICH  INCLUDES  DUMMIES*  * 

2  33H  EXCEEDS  THE  COMMON  STORAGE  BLOCK  *15) 

STOP 

780  CONTINUE 

C  *****  PRINT  THE  NGTOG  ARRAY 
DO  1100  1=1 *NOG 

WRITE! 6*3000) I *NAMEG( I ) *NGTO( I ) * (NGTOGd  * J) * J=1 *  10) 

3000  FORMAT! 1417) 

1100  CONTINUE 
RETURN 
END 
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Subroutine  INTRPL 


Subroutine  INTRPL  as  shown  in  Figure  1 1-4  is  called  from  INDATA  only  if 
there  are  NYC  reservoirs.  It  is  called  twice,  the  first  time  to  define  the 
storage  level  dividing  line  between  Normal  and  Drought  Warning  conditions,  the 
second  time  to  define  the  storage  level  dividing  line  between  Drought  Warning 
and  Drought  conditions. 

Each  time  it  is  called,  it  reads  the  beginning,  ending  and  inflection 
points  of  the  line  between  each  set  of  conditions.  It  then  determines  for  each 
day  of  the  year  the  storage  level  that  changes  from  Normal  to  Drought  Warning 
and  from  Drought  Warning  to  Drought  and,  of  course,  visa  versa.  INTRPL  returns 
to  INDATA. 
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Figure  1 1 -4.  Flowchart  for  Subroutine  INTRPL 


o  o  o  o 


SUBROUTINE  INTRPL(X) 

c  mmmmmmmmmmmmtmsmmmmmmmmmm 

THIS  SUBROUTINE  FIXES  THE  STORAGE  LEVELS  FOR  DETERMINING  THE 
SYSTEM  CONDITIONSt  USED  TO  OPERATE  THE  NYC  RESERVOIRS  AND  THE 
D  l  R  CANAL. 

mmmmmmmmmmmtmmttmmmmmmtttmsm* 

DIMENSION  X(31 f 12) . FLOW < 10) rNDAY ( 10) fMOKNT (13) 

DATA  MOKNT  /0»31 » 29.31 r 30. 31 ,20, 31 r 31 » 30»31 »30»31/ 

IM011=1 
IM01=2 
IDY1=1 

READ (5  > 100 ) FC » NCARDS 
IF<FC»NE. 1HF . AND.FC.NE. 1HG>  GO  TO  2100 
MR1TE<6» 101 )FC» NCARDS 
100  F0RMAT<A1f4X»I10> 

101  FORMAT (lX'Ali4XrI10) 

READ ( 5 ? 200 )  NMOiNDY ,FLOW< 1 ) 

URITE<6»200)  NMO»NDY»FLOW( 1 ) 

FLOW( 1 )=FLOW( 1 XI 550. 

X(1 >l)=FLOU(l) 

NDAY ( 1 )=MGKNT(NMO)+NDY 
200  FORMAT (5Xi2I5fF10.0) 

DO  1000  I=2»NCARDS 
READ(5i200)  NMO > NDY  fFL0W( I ) 

WRITE(6»200)  NMOfNDYtFLOM(I) 

FLOW ( I ) =FLOW( I ) *1550 . 

M0NTHS=0 
DO  300  II=1»NM0 
MONTHS=MONTHS+MOKNT ( I I ) 

300  CONTINUE 

NDAY ( I ) =NONTHS+NDY 
11=1-1 

DELDAY=NDAY(I)-NDAY(II) 

DELFLO=FLOW<I)-FLOW<II) 

SLOPE=DELFLO/DELDAY 
I IDY=DELDAY 
DO  500  IDY=1 f I IDY 
IDY1=IDY1+1 

IF ( IDY1 .LE.MOKNT ( IM01 ) )G0  TO  400 

IM01=IMOl+l 

IDY1=1 

IFlIMOl »6T • 13)IM01=2 
IMO 1 1=1 M01 — 1 
400  CONTINUE 

X< IDY1 » IM01 1 >=IDY*SLOPE+FLOW< I I > 

500  CONTINUE 
1000  CONTINUE 
RETURN 

2100  URITE(6t 2101 ) 

2101  FORMAT (57H  ERROR  1  INCORRECT  *Fa  OR  "G*  CARD  !  PROGRAM  TERMINATED 
.  ) 

STOP 

END 
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Subroutine  LINKER 


Subroutine  LINKER  is  called  from  R6FL  to  link  the  model  nodes  to  the 
incremental  inflow  tape  ( I  IT)  nodes.  Its  flow  chart  is  given  in  Figure  1 1 -5 . 

It  is  essential  that  the  model  node  numbers  correspond  to  the  tape  node  numbers 
that  are  given  in  Table  1 1 1-2 .  The  first  block  of  information  on  the  incremen¬ 
tal  inflow  tape  contains  the  node  numbers.  These  are  compared  with  all  of  the 
specified  node  numbers  in  the  input  data.  If  any  do  not  match,  a  warning  is 
printed,  containing  the  node  number  of  the  unclaimed  incremental  inflow,  on  the 
tape. 


Some  model  nodes  are  special  and  need  to  be  flagged.  They  are  Merrll 
Creek,  a  ficticious  Belvidere  where  Merrill  Creek  skims  and  the  D  &  R  canal. 
Each  of  these  has  a  special,  hard-coded  operation  performed  on  it  and  therefore 
needs  to  be  flagged.  LINKER  returns  to  RGFL. 


Figure  1 1 -5.  Flowchart  for  Subroutine  LINKER. 

(Continued) 
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SUBROUTINE  LINKER 

COHHON  /CONTROL/  NYEAR.NTAPE.NOUT .NOG.NNQDE. ITAPE. IYR. INO.IDY 
COMMON  /CHANL/  NGTO( 100) . NAHEG< 100) .NGTOGI 100. 10) . NOUTPT(IOO) » 

.  LINK ( 100) *Q( 100) *QIN( 100) *QFK(4) >QFK0(4) >QPR(4) . 

.  QPRQ<4) 

COMMON  /FLAGS/  NTR.NDR.MER.MPR.MTG.MRES.NBEL.MBL 
DIMENSION  NODE(IOO) 

c  mtmmmmmmmmmummmmmmmmummmm 

THIS  SUBROUTINE  SORTS  THE  INFLOW  INFORMATION*  ASSIGNS  IT  TO  THE 
APPROPRIATE  CHANNEL  AND  PRESERVES  THE  ORDER  FOR  REFERENCE 
IN  SUBROUTINE  INFLOW, 

NTR=0 
NDR=0 
MTG=0 
MER=0 
NBEL=0 

READ ( ITAPE .50) NNODE 
READ (I  TAPE *50) ( NODEIN ) *N=1» NNODE) 

50  FORMAT (15) 

DO  100  N=1 *NNOBE 
DO  200  NN=1 *  NOG 

IF  (NODE(N) . EQ.NAMEG(NN) )G0  TO  90 
IF (NAMEG(NN) ,EQ, 1155)  MER=NN 
IF ( NAMEGCNN ) . EG • 1 151 )  NBEL=NN 
200  CONTINUE 

WRITE(6*2000)N0DE(N) 

2000  FORMAT ( 30H  WARNING i  NO  MATCH  FOR  NODE  *14) 

GO  TO  100 
90  LINK(NN)=N 

IF<N0DE(N).EQ.1310)NTR=NN 
IF ( NODE (N) ,EQ, 1300)NDRSNN 
IF (NODE<N) ,EQ, 1 130 )MTG=NN 
100  CONTINUE 

WRITE (6* 150)  NTR,NDR.NER*HPR*NTG*NBEL.MBL 
150  FORMAT (///* 1X*27HTHE  TRENTON  SLOT  NUMBER  IS  *14* 

.  /.1X.29HTHE  DIR  CANAL  SLOT  NUMBER  IS  *14* 

.  /*1X*29HTHE  MERRILL  CR  SLOT  NUMBER  IS. 14* 

.  /.1X.29HTHE  PROMPTON  SLOT  NUMBER  IS  .14. 

.  /» 1X.29HTHE  MONTAGUE  SLOT  NUMBER  IS  .14. 

.  /.1X.29HTHE  BELVIDERE  SLOT  NUMBER  IS  .14. 

,  /. 1X.29HTHE  BLUE  MARSH  SLOT  NUMBER  IS. 14) 

RETURN 

END 


11-29 


Subroutine  INFLOW 


Subroutine  INFLOW  is  the  first  subroutine  called  within  the  daily  loop. 
Its  purpose  is  to  link  today's  incremental  inflows  to  the  proper  model  nodes. 
Because  the  incremental  inflows  are  input  in  the  US6S  format  (station  number 
followed  by  eight  daily  flows),  a  special  routine  is  required.  Figure  I I -6 
shows  the  basic  flowchart  for  INFLOW. 

An  8-day  counter  is  checked  to  see  whether  the  next  8  days  of  inflows 
should  be  read.  Once  read  in,  today’s  inflows  are  set  and  linked  to  the  proper 

node.  If  the  counter  is  not  at  8,  there  are  three  possibilities.  One,  if  it  is 

during  the  midweek  then  simply  today's  inflows  are  set.  Two,  if  it  is  the  last 
day  of  the  month,  the  daily,  monthly  and  yearly  counters  must  be  incremented 
properly  and  then  today's  inflows  are  set.  Three,  if  the  model  is  in  backup, 
the  counter  will  not  be  8.  However,  backup  is  not  checked  unless  it  is  the  last 

day  of  the  month.  If  both  backup  is  occurring  and  it  is  the  last  day  of  the 

month,  the  next  8  days  of  inflows  are  read  in  without  incrementing  the  daily, 
monthly  and  yearly  counters. 

It  is  possible,  during  a  50-year  simulation  to  come  to  the  end  of  the 
Incremental  Inflow  Tape  ( I  IT)  before  the  daily  and  yearly  counters  are  finished. 
If  this  is  the  case,  the  simulation  is  ended  with  the  message  -  END  OF  INPUT 
ENCOUNTERED  IN  INFLOW  -  and  the  output  subroutines  are  called.  Otherwise, 

INFLOW  returns  to  RGFL. 
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Figure  II-6.  Flowchart  for  Subroutine 
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Figure  1 1-6.  Flowchart  for  Subroutine  INFLOW. 

(Continued) 
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SUBROUTINE  INFLOU 

********************************************************************** 

INFLOW  PROVIDES  DAILY  INCREMENTAL  INFLOW  TO  ALL  MODEL  NODES 
AND  OPERATES  WALLENPAUPACK  WHEN  SPECIFIED  TO  DO  SO 
*********************************************************************** 
COMMON  /CONTROL/  NYEAR»NTAPF»N0UT»N06» NNODE»ITAPE»IYRtIMOi IDY 
COMMON  /CHANL/  NGTO( 100) .NAMEGUOO) »N6T0G( 100»10> fNOUTPTI 100) » 

.  LINK ( 100) fQ< 100) f OIN( 100) rDFK(4) rQFK0(4) »QPR(4) » 

.  QPR0(4) 

COMMON  /RESERV/  NRES>XFLOW» IFLA6»V0L( 10) .LAGRES(IO) tLNKRES(lO) » 

.  EVAP< 12*10) .STORE! 10) »CONS( 12» 10) *Q01 ( 100) * 

.  002(100) »NOGG* ICKrFLOWLrZFLOW? ICK1 

COMMON  /CAL/  CALIB1 ,JOUT» JMIN 
COMMON  /LAG/  QOUT( 100 ) »XM( 100) » B< 100) 

COMMON  /STOR/  ST0RE1 ( 10) »ST0RE2< 10) »ST0RE3( 10) »ST0RE4( 10) 

COMMON  /DATE/  I YR3f IM03» IDY3» I YR2r IM02» IDY2» I YR1 » IM01 » IDY1 » 

.  Q1 ( 100 ) > Q2( 100) 

COMMON  /PADER/  ACONS( 12» 10) »TRENT » YFLOW»REL( 10) »NSTEP»MDATE 
COMMON  /BACKUP/  YF5» YF4» YF3» YF2»YF1 
COMMON  /WALLIN/  MWL 
DIMENSION  WALL (60) 

DIMENSION  FLOW( 10i65) »KNT (12) 

C  *****  WALL  CONTAINS  THE  STEPPED  UP  RELEASES  FROM  WALLENPAUPACK 
C  DURING  THE  1960'S  DROUGHT 

DATA  WALL/39. 4»355. »291. ,424. » 549. » 741. »70.»70.»70.»70.»  0.il8.4» 

1  365. *404. >243.  »407. »346.  »299. »70. »70. »70. »70. »415. »484. » 

2  302. »382. >236*  >384. >470. »19B. »70. »70. »70. »70. *  0.»  0.» 

3  0. » 145. »320. »460. »494.»313.»70.»70.»70.»70.»  0.»8.36» 

4  390. *357. ,277. ,433. , 422. ,570. , 70. ,70. ,70. ,70. ,357. ,564./ 
DATA  KNT/7»4»7»6»7»6»7>7,6,7,6»7/ 

C  *****  INCREMENT  THE  DAY  OF  THE  WEEK  BEING  PROCESSED 
KOUNT=KGUNT+ 1 

C  *****  IF  KOUNT  =  8  A  NEW  WEEKS  CARDS  MUST  BE  READ 
IF (KOUNT-8 ) 4 ,5 ,5 

4  CONTINUE 
IDAY=IDY-24 

C  *****  CHECK  FOR  THE  LAST  DAY  OF  THE  MONTH  AND  FOR  BACKUP  OCCURRANCE 
IF( IDAY-KNT  ( IMO) ) 30, 1 , 3 

1  IFUMO-2)  3,2,3 

2  IF(MQD( IYRi 4 ) .EQ.O)GO  TO  30 

3  IF ( IFLAG)5»  21  * 6 

6  IF(IFLAG-l)  5,21,5 

5  CONTINUE 

C  *****  READ  INFLOWS  FOR  THE  NEXT  WEEK  AND  SAVE  THE  INFLOWS  OF  THE 
C  *****  LAST  2  DAYS  OF  THE  PREVIOUS  WEEK 
DO  10  J=1 *NNODE 
K9*K0UNT+1 
K10=K0UNT+2 

C  *****  SAVE  THE  LAST  2  DAYS  VALUES  IN  THE  FLOW  ARRAY 
FLOW( 1 , J)*FL0W(K9, J) 

FL0W(2» J)»FL0W(K10> J) 

10  CONTINUE 

11  CONTINUE 

C  *****  READ  ANOTHER  WEEKS  FLOWS 
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DO  20  J-l f NNODE 

READ ( I TAPE* 15) IYR. IMO  » IDY.  (FLGWd.  J) i 1*3. 10) 

IF (EOF ( I T APE )) 2000 >20 
15  FORMAT (18X.3I2.8F7.0) 

20  CONTINUE 

IF (MWL)23* 23.223 

C  *****  FIX  TO  REPLACE  WALLENPAUPACK  RELEASES  FROM  JUNE'62  THROUGH  MAY'67 
223  NDATE=IYR*100+IN0 

IF (NDATE.LT. 6206. OR. NDATE.GT. 6705)  60  TO  23 
C  *****  INCREMENT  THE  MONTHLY  REPLACEMENTS 
IF(NDATE.NE.MDATE)  NSTEP=NSTEPf 1 
DO  201  I=3f 10 
FLOU( I » 1 1 >=WALL(NSTEP) 

201  CONTINUE 
MDATE=NDATE 
GO  TO  23 

21  CONTINUE 

C  *****  INCREMENT  THE  MONTHS  DURING  THE  BACKUP  PROCEDURE 
IM0=IN0fl 
IDY=1 

IF(IM0-13)23. 22.23 

22  CONTINUE 
IM0=1 

I YR=IYR+1 

23  CONTINUE 

C  *****  RESET  THE  COUNTER  AFTER  READING  A  NEU  WEEKS  FLOWS 
K0UNT=0 
IDY= ( IDY-1 )*8 
GO  TO  30 

C  *****  ENTRY  POINT  FOR  BACKUP  PROCEDURE 
ENTRY  NEWFLO 
IF (NRES.EQ.O)RETURN 
WRITE! 6 .3000) IYR* IMO. IDY 
3000  FORMAT (22H  LOW  FLOW  AT  MONTAGUE  .312) 

C  *****  RESET  THE  DATE  TO  THREE  DAYS  AGO 
IDY=IDY3-1 
IDY1=IDY3 
IM01=1M03 
IYR1=IYR3 
I M0=IMO3 
I YR= I YR3 
K0UNT=K0UNT-3 

C  *****  RESET  TRENTON  TARGET  TO  THREE  DAYS  AGO 
YFL0W=YF3 
YF l  =  YF4 
YF2=YF5 
YF3=YF6 

C  *****  CALCULATE  XFLOW— THE  RELEASE  TO  MONTAGUE 
XFL0W*CALIB1*(FL0WL-Q( ICK ) ) 

C  *****  SET  IFLAG  TO  COUNT  THE  DAYS  UNTIL  THE  BACKUP  IS  OVER 
IFLAG*3 

C  *****  RESET  THE  FLOWS  TO  THREE  DAYS  AGO 
DO  25  KS1 .NOGG 
01 (K)*Q02(K) 
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GOUT  <K)=002(K) 

25  CONTINUE 

C  *****  RESET  THE  RESERVOIR  LEVELS  TO  THREE  DAYS  AGO 
DO  27  K=1*NRES 
ST0RE1(K)*ST0RE4(K) 

STORE (K)=ST0RE4(K) 

27  CONTINUE 

C  *****  SET  BLUE  HARSH  HAXIHUH  VOLUME  TO  THE  CORRECT  LEVEI 
C  IFUN0.EQ.9)  VOL(NRES>=11545. 

C  IHIM0.EQ.3)  VOL <NRES> =8873. 

RETURN 
30  CONTINUE 

C  *****  INCREMENT  FOR  THE  NEXT  DAYS  INFLOW 
IDY=I0Y+1 
IF'0INT=K0UNT+3 
DO  40  J=1 *NNODE 
QIN< J)=FLOW( IP01NT *  J) 

40  CONTINUE 

1001  CONTINUE 

C  *****  LINK  THE  NEW  INFLOWS  TO  THE  PROPER  NODE 
DO  100  1=1 »NOOG 
0(I)=0. 

N=LINK( I ) 

Q( I ) =OIN(N) 

100  CONTINUE 
RETURN 

2000  WRITER. 2100) 

2100  F0RHAT(35H  END  OF  INPUT  ENCOUNTERED  IN  INFLOW  ) 

IF (JOUT >2200*2300 *2400 
2200  CALL  OUTPUT 
GO  TO  2500 
2300  CALL  MONTH 
GO  TO  2500 
2400  CALL  HONTH 
CALL  OUTPUT 
2500  CONTINUE 
STOP 
END 
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Subroutine  ROUTE 


Before  all  of  the  flows  can  be  routed,  the  reservoir  releases  must  be 
calculated.  As  shown  in  Figure  1 1-7,  ROUTE  calls  RESOP  immediately  after 
defining  the  lag  functions.  Next,  any  consumptive  losses  are  subtracted  from 
the  appropriate  node.  Finally,  the  flows  are  lagged  and  added  from  upstream  to 
downstream  and  at  each  downstream  node,  the  incremental  inflows  are  added  in. 
ROUTE  returns  to  RGFL.  Immediately  upon  return,  the  flows  are  written  on  an 
output  tape,  providing  the  model  is  not  backing  up. 


11-36 


r>  o  o  r>  r> 


SUBROUTINE  ROUTE 

c  *********************************************************************** 

SUBROUTINE  ROUTE  LAGS  THE  FLOWS  FROM  UPSTREAM  TO  DOWNSTREAM 
AND  ADDS  IN  THE  PROPER  INCREMENTAL  INFLOWS.  ROUTE  CALLS  RESOP 
TO  OPERATE  THE  RESERVOIRS  IN  ORDER  TO  DETERMINE  THEIR  DOWNSTREAM 


FLOWS 

*********************************************************************** 

COMMON  /CONTROL/  NYEAR.NTAPE.NOUT  .NQG.NNODE. ITAPE. IYR. IMO. IDY 
COMMON  /CHANL/  NGTOt 100) »NAMEG< 100) »NGTOG< 100. 10) »NOUTPT< 100) » 

.  LINK(100)»Q<100>,QIN<100).QFK<4),QFK0(4)»QPR<4>. 

.  QPR0(4) 


COMMON  /RESERV/  NRES » XFLOW > IFLAG » VOL (10)* LA6RES ( 10 ) » LNKRES ( 10 ) » 
EVAP( 12i 10) t STORE ( 10) .CONS( 12. 10) .Q01 ( 100) r 
Q02< 100) .NOGG » ICK.FLOWL. ZFLOW. ICK1 
COMMON  /LAG/  QOUT< 100) .XM< 100) »B< 100) 


COMMON  /CONSUM/  NCON. IC0N<50) »CLOSS( 10. 12) 
C  *****  DEFINE  THE  LAG  FUNCTIONS 


X  LAG  < I ) =Q  < I )*XM( I )+B< I ) 

YLAG( I ) =QOUT (I)*XM<I)+B(I) 

C  *****  CALL  RESOF - 0FERATE  THE  RESERVOIRS 

CALL  RESOF 

C  *****  SUBTRACT  THE  CONSUMPTIVE  LOSSES  FROM  SELECTED  NODES 
IF (NCON.EQ.O)  GO  TO  951 
DO  950  1  =  1 » NCON 
I I=ICON( I ) 

0 ( 1 1 ) =Q< 1 1 ) -CLOSS< I . IMO ) 


950  CONTINUE 

951  CONTINUE 

C  *****  FIND  THE  UPSTREAM  NODES  FOR  EACH  NODE  AND  LAG  AND 
C  *****  ADO  THE  FLOWS  FROM  THE  UPSTREAM  NODES 
DO  900  1=1 . NOGG 
DO  100  11=1.10 

IF (NGTOG( 1 . 1 1 ) .EQ< 0 )G0  TO  101 
1 1 I=NGTOG< I .II) 

Q( I )=Q  < I ) +QOUT (III )*YLAG< III)+Q(III)*(1.0-XLAG(III)) 

100  CONTINUE 

101  CONTINUE 
900  CONTINUE 

DO  1000  1*1. NOGG 

C  *****  0(1)  IS  THE  DAILY  FLOW  AT  EACH  NODE 

C  *****  QOUT(I)  SAVES  THE  DAILY  FLOW  ARRAY  FOR  THE  NEXT  DAYS  PROCESSING 
OOUT ( I )=Q( I ) 

1000  CONTINUE 
RETURN 
END 


Subroutine  RESOP 


Subroutine  RESOP  whose  flowchart  is  given  in  Figure  1 1 -8  operates  the 
reservoirs  and  makes  releases  to  augment  both  the  Montague  and  Trenton  flows. 

Upon  entering  RESOP,  the  beginning  of  the  month  is  checked.  If  it  is 
the  first  day  of  the  month,  three  things  occur:  1)  If  Blue  Marsh  is  included, 
its  maximum  storage  is  changed  because  a  higher  reservoir  level  is  maintained 
during  the  summer  months  for  recreation  purposes;  2)  The  Special  Outputs  that 
require  the  beginning  storage  levels  are  set  for  all  reservoirs  to  the  current 
storage  and;  3)  The  previous  month's  minimum  storage  is  saved  for  output  pur¬ 
poses. 


On  every  new  day,  the  System  Condition  is  determined  using  the  total 
beginning  NYC  reservoir  storage.  The  condition  will  not  be  changed  during 
backup.  The  System  Condition  then  specifies  the  amount  of  the  diversions  and 
the  objectives. 

The  next  logical  step  is  to  calculate  new  storages  using  inflow,  con¬ 
servation  releases  and  evaporation  factors.  Before  this  can  be  done  the  NYC 
evaporation  needs  to  be  set.  8ecause  the  net  inflow  values  are  available  for 

the  NYC  reservoirs  after  they  are  put  on  line,  evaporation  must  be  set  to  zero 

for  that  particular  reservoir  at  the  appropriate  date.  Once  this  is  done,  the 
new  storages  for  all  the  reservoirs  are  calculated. 

First  the  NYC  diversion  is  taken  from  the  fullest  NYC  reservoir  and  its 
new  storage  calculated.  Then  any  release  needed  to  maintain  the  Montague  flow 
is  taken  from  the  fullest  NYC  reservoir,  using  the  new  storage  just  calculated. 
Next,  any  release  needed  to  augment  Trenton  is  taken  from  the  fullest  lower 
basin  reservoir. 

Merrill  Creek  is  an  optional  reservoir.  If  it  is  online,  the  Trenton 

flow  is  checked:  if  the  flow  is  less  than  3100  cfs,  Merrill  Creek  releases 

126.9  cfs  to  Riegelsvil le.  If  the  flow  is  greater  than  3200  cfs,  Merrill  Creek 
skims  from  Belvidere  to  replenish  its  storage. 
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Prompton  reservoir  being  located  above  Montague  will  change  the  NYC 
reservoirs'  releases  if  the  Montague  objective  does  not  also  change.  At  this 
point  in  RESOP,  this  is  done.  The  natural  inflow  to  Prompton  is  routed  down  to 
Montague  without  taking  into  account  any  incremental  inflows  from  further  down 
the  basin  using  the  actual  route  times.  The  total  releases  from  Prompton  reser¬ 
voir  are  then  also  routed  in  the  same  manner.  The  difference  is  determined  and 
then  added  to  the  Montague  target.  A  positive  change  indicates  that  the 
Prompton  reservoir  releases  are  greater  than  the  natural  Prompton  flow.  The 
objective  increases  so  that  these  releases  move  on  down  to  Trenton  and  are  not 
used  to  augment  Montague's  flow.  A  negative  change  indicates  the  opposite: 
Prompton' s  natural  flow  was  greater  than  the  reservoir  releases  and  therefore 
the  model  reduces  the  Montague  target  objective  to  keep  the  reservoirs  from 
releasing  more  than  they  should  under  conditions  when  Prompton  was  not  included. 

The  final  function  of  RESOP  is  to  save  the  Special  Outputs.  These 
Special  Outputs  are  different  from  the  flows  at  any  particular  node.  They  con¬ 
sist  of  total  storages,  releases  neded  to  maintain  Montague  or  Trenton,  minimum 
monthly  storages  and  other  various  requests.  They  are  hard-coded  into  the 
program  but  can  be  changed  easily.  RESOP  returns  to  ROUTE. 


Figure  II-8.  Flowchart  for  Subroutine  RESOP 
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Figure  1 1-8.  Flowchart  for  Subroutine  RESOP. 

(Continued) 
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Figure  1 1  - 8 .  Flowchart  for  Subroutine  RESOP. 

(Continued) 


Figure  1 1 -8 .  Flowchart  for  Subroutine  RESOP. 

(Continued) 
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Figure  1 1 -8.  Flowchart  for  Subroutine  RESOP 

(Continued) 
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SUBROUTINE  RESOP 

C  RESOP  OPERATES  THE  RESERVOIRS  IN  THE  MODEL 

C  1)  CONSERVATION  RELEASES  ARE  MADE  — BASIC  OR  AUGMENTED 

C  2)  RELEASES  ARE  MADE  AND  DIVERSIONS  AND  EVAPORATION 

C  ARE  SUBTRACTED  EACH  DAY 

c  *********************************************************************** 

COMMON  /CONTROL/  NYEAR . NTAPE.NOUT .NOG.NNODE. ITAPE . I YR» IMO. IDY 
COMMON  /RESERV/  NRES. XFLOW . IFLAG. VOL ( 10 > .LAGRES ( 10) iLNKRESt 10) » 

.  EVAP( 12i 10) » STORE! 10)  »CONS( 12. 10) .001 ( 100) » 

.  002(100) .N0GG.ICK.FL0WL.ZFL0W.ICK1 

COMMON  /CHANL/  NGTO( 100) »NAMEG( 100) »NGTOG( 100.10) .NOUTPT ( 100) » 

.  LINK ( 100) *G( 100) >QIN( 100) »QFK(4) . QFK0(4 ) »QPR(4) » 

.  QPR0(4) 

COMMON  /PHAS2/  NYCRES . ICHK»STOREN( 31 . 12) .STORED (31 >12) t 
.  DI VERS  <  3 ) f OBJECT ( 3 ) .  EVP( 12.10). I EVP  *  DFLOM( 3 )  . 

.  TARG<3) 

COMMON  /PADER/  ACONS < 12 » 10) . TRENT » YFLOM» REL ( 10) » NSTEPi MDATE 

COMMON  /FLAGS/  NTR»NDR»MERcMPR»MTGf MRESr NBELrMBL 

COMMON  /DATE/  I YR3 . IM03 » IDY3 » I YR2 » IM02. IDY2* I YR1 r IM01 » ID Y1 » 

.  QK100)  .02(100) 

COMMON  /MOMIN/  QMM< 12.50. 10) »QMIN< 10) 

DIMENSION  WIT(10).IDATE(10) .QBS(10> 

DIMENSION  FXM(3).FB(3).SPL(10) 

DATA  FXM/-1 . 4221E-05.-6 .7208E-06 .-4 .8167E-06/ 

DATA  FB/. 0833. .2083.. 4583/ 

C  *****  IDATE  CONTAINS  THE  SWITCH  DATES  FOR  TURNING  ON  AND  OFF 
C  NYC  RESERVOIR  EVAPORATION  WITH  THE  NYC  INFLOW  DATA 

DATA  I DATE/530701. 540101 .540201 .540917.550601 .550701 . 

.  590801 » 600101.631001 .771001/ 

C  *****  SET  LAG  FUNCTIONS  FOR  DTARG  CALCULATIONS 
FKIX< I ) =QFK ( I )*FXM(I )+FB(I ) 

FKLR(I)=QPR(I)*FXM(I)+FB(I) 

FKl Y( I )=QFKO( I >*FXM( I )+FB( I ) 

FKRY( I ) =QPRO( I )*FXM( I ) fFB( I ) 

C  *****  IF  THERE  ARE  NO  RESERVOIRS.  RETURN  TO  MAIN  PROGRAM 
IF  (NRES.EQ.O)  RETURN 
NRESlrNRES 

C  *****  FOR  BLUE  MARCH  INCLUSION.  MUST  NOT  USE  IT  TO  MEET  TRENTON  FLOW 
IF (MBL)2»2» 1 
1  NRES1=NRL3-1 

C  *****  IF  IN  BACKUP.  SKIP  THE  RESETTING  OF  QMIN  AND  BLUE  MARSH 
IF(IFLAG.GT.O)  GO  TO  10 

*****  FIND  THE  BEGINNING  MONTH  STORAGE  AND  START  SEARCHING  FOR  MINIMUM 
IF  < IDY  —  1 )  10.20.10 

*****  CHANGE  BLUE  MARSH  STORAGE  LEVEL  FROM  SUMMER  TO  WINTER 
0  IF ( MBL  >45.45.21 

1  IF(IM0.LT.4.0R.IM0.GT.9)  GO  TO  40 

VOL ( NRES  >  =  1 1545 . 

GO  TO  45 

40  VOL ( NRES ) =8873 . 

45  DO  25  1=1. NRES 

OBS< I >=STORE( I ) 

IF ( I MO- 11 >50.55.55 
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IMF=IM0+2 
IYF=I YR-27 
GO  TO  35 
55  IMF=IM0-10 

I YF=IYR-26 

35  IF (IMF ) 10> 10*36 

36  QMM < IMF  *  I YF  *  I >  =QM IN < I ) *0 . 00064627 

25  QMIN(I >=9999999 . 

10  CONTINUE 

DO  30  I=1*NRES 
SPL ( I ) =0 . 

30  QMIN< I )=AMIN1 (QMIN( I ) » STORE ( I )  > 

XNYC=0. 

QT0T=0. 

GV0L=0. 

QTOT  A=0 « 

F'RIN=Q<MPR) 

lF(NYCRES.LT.l)  GO  TO  101 

C  mil  COMPUTE  TOTAL  STORAGE  AVAILABLE  IN  THE  NYC  RESERVOIRS 
DO  100  N=1 f NYCRES 
100  QTOT=OTOHSTORE(N) 

101  CONTINUE 

IF(NRES.EQ. NYCRES)  GO  TO  104 

C  *****  COMPUTE  TOTAL  STORAGE  AVAILABLE  FOR  ALL  OTHER  RESERVOIRS 
KRES=NYCRESH 
IF(MER.NE.O)  KRES=KRES+1 
DO  103  N=KRES*NRES1 
QTOTA=QTOTA+STORE(N) 

QVOL=VOL ( N ) +QVOL 

103  CONTINUE 

104  CONTINUE 

IF ( IFLAG.GE .0)G0  TO  205 

IOLD=ICHK 

ICHK =2 

IF ( NYCRES .LT . 1 )  GO  TO  1001 
IF(QTOT .GT . STOREN< IDY> IMO) )  ICHK=1 
IF ( QTOT .LT . STORED ( I DY » INO) )  ICHK»3 
205  CONTINUE 

FLOWL=OBJECT ( ICHK ) 

MITHH=DI VERS( ICHK) 

TRENT=TARG( ICHK) 

Q(NDR)=-DFLON( ICHK) 

C  *****  SKIP  SWITCHING  EVAPORATION  ON  AND  OFF  IF  NYCRES  DOES  NOT  EOUAL  3 
IF(NYCRES.NE.3>  GO  TO  420 
IF< IYR.LT.53.0R. IYR.GT.63)  GO  TO  420 
NDATEsIYR*10000+!M0*100+IDY 

IF (NDATE .EQ. I  DATE ( I EVP) )  GO  TO  (310» 330»310»350»370»350» 

.  330 » 3 10 » 390 » 420) r I EVP 

GO  TO  420 

C  *****  EVAPORATION  SWITCH  FOR  NEVERSINK 
310  DO  320  IMS1 r 12 
EVAP<IM»3)=0. 

EVAP( IM» 1 )*0. 

320  CONTINUE 
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330  DO  340  IH=1 *  12 

EVAP< IM»3)=EVP< I M » 3 ) 

340  CONTINUE 
GO  TO  410 

C  *****  EVAPORATION  SWITCH  FOR  PEPACTON 
350  DO  360  I M= 1 *12 
EVAP<IM»1)=0. 

360  CONTINUE 
GO  TO  410 
370  DO  380  IM=1 . 12 

EVAP(IM.l)=EVP<IM.l) 

380  CONTINUE 
GO  TO  410 

C  *****  EVAPORATION  SWITCH  FOR  CANNONSVILLE 
390  DO  400  IM=1 » 12 
EVAP( IM»2)=0. 

400  CONTINUE 

410  CONTINUE 

IF ( IFLAG.GE.O)  GO  TO  420 
WR ITE  <  6 » 41 1 )  IEVP.NDATE 

411  FORMAT (24H  EVAPORATION  SWITCH  NO.  r 1 2 r 1 1H  THROWN  ON  .16) 
IEVP=I£VP+1 

420  CONTINUE 

DO  430  1=1 iNRES 
430  WIT< I )=0. 

QSAV=-300. 

DO  1000  11=1 .NYCRES 
IX=NOGG+II 

IF(Q( IX) .LT. QSAV)  GO  TO  1000 
QSAV=Q( IX) 

NSAV=I I 
1000  CONTINUE 
1001  CONTINUE 

UIT (NSAV)=WITHH 
DO  2000  N=1 *NRES 
NN=LNKRES<N) 

STORE(N)=STORE<N)FQ<NN) 

STORE <N)=STORE(N)-EVAP(INO»N)-WIT(N> 

2000  CONTINUE 

IF(ICHK.EQ.l)  GO  TO  3001 
C  *****  MAKE  BASIC  CONSERVATION  RELEASES 
DO  3000  N=1 .NRES 
NN=LNKRES(N) 

IN=NOGG+N 

Q(NN)=STORE<N)-VOL<N) 

IF (Q(NN)-CONS( IMO.N) )  2100.2400.2300 
2100  Q(NN)=CONS( IMO.N) 

GO  TO  2400 

C  *****  OPTIONAL  OUTPUT  \  SPILLS  ARE  COMPUTED  FOR  EVERY  RESERVOIR 
300  SPL(N)=Q(NN)-CONS< IMO.N) 

400  STORE! N) “STORE (N)-Q(NN) 

Q(IN)=ST0RE(N)/V0L(N)*100. 

3000  CONTINUE 


N)  rj 


80  TO  121 

3001  CONTINUE 

C  *****  HAKE  AUGMENTED  RELEASES 
DO  3002  N=1 *NRES 
NN=LNKRES(N) 

IN=NQGG+N 

n(NN)=STORE(N)-VOL(N> 

IF (Q ( NN ) -ACONS< IHO.N) )  2500 . 2800 t 2700 
2500  G(NN)=ACONS( INO.N) 

GO  TO  2800 

700  SPL ( N ) =Q  <  NN ) -ACONS (IHO.N) 

800  STORE(N)=STOR£(N)-Q(NN) 

G ( IN) “STORE (N)/V0L(N>*100. 

3002  CONTINUE 
121  CONTINUE 

DO  3004  1  =  1 .  NRES 

3004  REL( I )=0. 

C  *****  BEGIN  CALCULATION  OF  RELEASE  REQUIRED  TO  MAINTAIN 
C  *****  2ND  TARGET'S  FLOW 

IF < YFLOM )  3008.3008*3005 

3005  CONTINUE 
0SAV=-300. 

DO  3006  N=KRES.NRES1 
NT=NOGG+N 

IF<G(NT) .LT.QSAV)  GO  TO  3006 

NFULL=N 

OSAV=G(NT) 

3006  CONTINUE 

REL(NFULL)=YFLOH 
NN=LAGRES( NFULL ) 

IN-N0G6+NFULL 

STORE  <  NFULL ) “STORE  <  NFULL ) -REL ( NFULL ) 

Q(NN)=Q(NN)+REL( NFULL) 

Q(IN)=STORE (NFULL) /VOL (NFULL >*100. 

3007  CONTINUE 

3008  CONTINUE 

C  *****  BEGIN  TO  COMPUTE  RELEASE  REQUIRED  TO  MAINTAIN  1ST  TARGET'S  FLOW 
IF(XFLOU >5000. 5000. 4100 
4100  QSAV=-300. 

DO  4000  11=1 .NYCRES 
IX=NOGG+II 

IF(Q(IX).LT.QSAV)GO  TO  4000 
QSAV=Q( IX  > 

ISAV=II 
4000  CONTINUE 

REL(ISAV)=XFLOM 
NN=LAGRES( ISAV) 

IN=NOGG+ISAV 
Q ( NN ) “XFL  OU+Q ( NN ) 

STORE! ISAV) “STORE! ISAV) -XFLOW 
Q( IN)=STORE( ISAV)/VOL( ISAV) *100. 

5000  CONTINUE 

C  *****  MERRILL  CREEK  OPTION 
IF (MER.EQ.O)  60  TO  5500 
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IF ( ZFLOW) 5510 t 5520 >5520 
5510  IN“NOGG+HRES 

IF(QdN).LE.O)  60  TO  5500 
NN“LNKRES(HR£S) 

Q<NN)=126.9+Q<NN) 

STORE (HRES) “STORE ( HRES )-Q(NN) 

Q ( IN ) =STORE  <  HRES ) /VOL ( HRES ) *100 . 

GO  TO  5500 
5520  IN=NOGG+HRES 

IF < 0 ( IN ) • GE < 100 • )  GO  TO  5500 
ZFLOW= ZFLOW- 100. 

IF (ZFL0W)5500>5500>5525 
5525  ZF1=V0L( HRES) -STORE (HRES) 

ZFL0W=AHIN1 ( ZFLOW >ZF1 ) 

ZFL0W=AHIN1 ( ZFLOW f 145 . ) 

STORE ( HRES ) “STORE ( HRES ) +ZFL0W 
Q( IN )=STORE< HRES) /VOL (HRES) >100. 
Q(NBEL)=Q(N»EL)-ZFLOW 
5500  CONTINUE 

C  *****  END  OF  HERRILL  CREEK  OPTION 
C  *****  BEGIN  PROHPTON  OPTION 

C  *****  FIRST  ROUTE  DOWN  THE  NATURAL  PROHPTON  FLOW 
IF (HPR) 5200 >5200 >5100 
5100  CONTINUE 

GFK( 1 )=PRIN 
DO  180  JK=1 >3 
K J=JKF1 

QFK ( K J ) =QFK  <  JK ) * ( 1 -FKLX  <  JK ) ) +QFKO ( JK ) *FKL Y ( JK ) 

180  QFKO ( JK ) “QFK < JK ) 

C  *****  NOW  ROUTE  THE  ACTUAL  RELEASES  FROH  PROHPTON 
QPR( 1 )=Q(HPR) 

DO  181  JK=1 >3 
KJ=JK+1 

QPR(K J)  =OPR(  JK )*(  1-FKLR(  JK)  )4QPR0(  JKXFKRY <  JK ) 

181  QPRO ( JK ) =QPR ( JK ) 

C  *****  SUBTRACT  NATURAL  FLOW  FROH  ACTUAL  RELEASES 
DTARG=QPR ( 4 ) -QFK ( 4 ) 

FLOWL=FLOWL+DTARG 
5200  CONTINUE 

C  *****  END  OF  PROHPTON  OPTION 

C  *****  BLUE  HARSH  OPTION  i  9  CFS  WATER  SUPPLY  WITHDRAWAL 
IF (NBL)5300>5300>5250 
5250  IN=NOGG+NRES 

IF(Q(IN).LE.O)  GO  TO  5300 
STORE (NRES) “STORE (NRES) -9.0 
Q (IN ) “STORE ( NRES ) /VOL  <  NRES ) 1 100 . 

5300  CONTINUE 

C  *****  END  BLUE  HARSH  OPERATION 
C  *****  SPECIAL  OUTPUTS 
NZ=N0GG+NRES+1 
C  *****  HONTAGUE  TARGET 
C  Q(NZ)=FLOWL 

C  NZ=NZ+1 

C  *****  RELEASES  TO  HAINTAIN  1ST  TARGET  NODE'S  FLOW 
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Q(HZ)*XFLOU 

NZ=NZ+1 

C  *****  RELEASES  TO  MAINTAIN  2ND  TARGET  NODE'S  FLOW 
Q(NZ)=YFLOW 
NZ=NZ+1 

C  *****  DIVERSIONS  TO  NYC  RESERVOIR  SYSTEM  IN  MGD 
Q(NZ)=WITHH/1 .5473 
NZ=NZ+1 

C  *****  NYC  RESERVOIRS'  TOTAL  STORAGE  IN  BILLION  GALLONS 
IF(NYCRES)  5320.5320.5310 
5310  Q(NZ)=QT0T*0. 00064627 

NZ=NZ+1 

C  *****  LOWER  BASIN  RESERVOIRS'  TOTAL  STORAGE SNOT  INCLUDING  MERRILL  CREEK 

C  *****  OR  BLUE  MARSH 

5320  Q(NZ>=QT0TA*0, 00064627 

C  *****  SPILLS  FROM  EACH  AND  ALL  RESERVOIRS 

C  DO  70  1=1 . NRES 

C  NZ=NZ+1 

C  Q(NZ)=SPL( I )-REL( I ) 

C70  IF(Q(NZ).LT.O.)  Q(NZ)=0. 

C  *****  BEGINNING  STORAGES  FOR  EACH  AND  ALL  RESERVOIRS 
C  DO  60  1=1, NRES 
C  NZ=NZ+1 

C60  Q(NZ)=QBS(I >*0.00064627 

C  *****  RELEASES  MADE  FROM  EACH  RESERVOIR  TO  MAINTAIN  TARGETS 
C  DO  4010  1=1, NRES 
C  NZ=NZ+1 
C  Q(NZ)=REL( I ) 

C4010  CONTINUE 

C  *****  NYC  DIVERSIONS  FROM  EACH  OF  THE  NYC  RESERVOIRS 
C  DO  4020  1=1,3 
C  NZ=NZ+1 
C  G(NZ)=WIT( I )/l .5473 
C4020  CONTINUE 
XFL0W=0 . 

ZFL0W=0 . 

RETURN 

END 


i 
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Subroutine  NEWFLO 


Subroutine  NEWFLO  is  called  from  RGFL  only  if  the  model  is  not  already 
in  backup  and  if  Montague  is  below  its  objective.  Figure  1 1-9  shows  the  simple 
structure  of  NEWFLO. 

NEWFLO  backs  the  simulation  up  three  days  and  reinitializes  the  stora¬ 
ges,  inflows  and  counters.  It  calculates  the  amount  of  release  needed  to  main¬ 
tain  the  objective.  It  then  multiplies  this  release  by  the  factor  used  to 
overcome  attenuation.  The  simulation  continues  as  before  only  no  flows  are 
written  on  tape  until  the  three-day  backup  is  completed.  The  model  is  set  up 
such  that  it  will  never  back  up  to  a  specific  day  more  than  once.  This  elimina¬ 
tes  the  possibility  of  releasing  twice  to  maintain  the  flow  on  any  particular 
day.  NEWFLO  returns  to  RGFL. 
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Subroutine  MONTH 


Figure  11-10  shows  the  primary  operation  of  Subroutine  MONTH.  The  data 
for  each  requested  output  must  be  linked  to  the  appropriate  node  or  Special 
Output.  If  an  output  node,  either  flows  and  Special  Outputs,  cannot  be  linked 
to  a  node  from  the  scratch  tape,  an  error  message  will  be  printed  ERROR  IN 
OUTPUT:  SELECTED  OUTPUT  NOT  FOUND.  This  will  result  in  a  table  of  zeros  but 
will  not  stop  the  printing  of  the  rest  of  the  tables.  Subroutine  MONTH  reads 
all  of  the  data  for  a  single  day  from  the  scratch  tape,  connects  the  proper  node 
to  output  request  and  begins  building  a  monthly  array.  Each  day,  the  values  are 
properly  connected  and  added  to  the  previous  sum.  At  the  end  of  each  month,  the 
monthly  average  is  calculated  and  stored. 

Once  all  of  the  data  has  been  read,  the  monthly  tables  are  printed. 

Each  node  has  its  coresponding  monthly  average  values  printed  by  year  in  rows 
and  by  month  in  columns.  The  overall  average  value  for  each  month  is  calcu¬ 
lated.  First,  the  flow  values  at  the  requested  nodes  are  printed.  Then  come 
the  Special  Outputs,  printed  in  the  order  specified  in  the  input  deck. 

When  the  printing  has  finished,  Subroutine  MONTH  prints  -  SIMULATION 
ENDED  NORMALLY  -  and  returns  to  RGFL. 
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SUBROUTINE  HONTH 

C  MONTHLY  OUTPUT  GENERATED  USING  DAILY  VALUES  AND  AVERAGES 

c  ********************************************************************** 

COMMON  /CONTROL/  NYEAR» NTAPE»NOUT»NOG» NNODE* ITAPE. IYRt IMO. IDY 
COMMON  /CHANL/  NGTO( 100) .NAMED ( 100) .NGTOG( 100. 10) .NOUTPT ( 100) f 

.  LINK ( 100) >Q( 100) rQIN( 100) »QFK( 4 ) »CFKO( 4 ) »QPR( 4 ) . 

.  QPRO( 4) 

COMMON  /INOUT/  NTITLE1 ( 10) »NTITLE2( 10) »NTITLE3( 10. 100) 

COMMON  /RESERV/  NRES.XFLOW. IFLAG»VOL( 10) »LAGRES( 10) »LNKRES( 10 > » 

.  EVAP< 12f 10) • STORE ( 10) .CONS( 12. 10) .001(100). 

.  Q02  ( 1 00 )  i  NOGG .  ICK » FLOML . ZFl.OW » ICK 1 

COMMON  /DATE/  IYR3. IM03. IDY3. I YR2. IM02* IDY2. IYR1 » IH01 » IDY1 » 

.  01(100) >02(100) 

COMMON  /MOMIN/  QMM( 12.50.10) 

COMMON  /CAL/  CALIB1 . JOUT, JMIN 

DIMENSION  FMAX( 12»75)»FMIN(12»75) »FMEAN( 12.75) » ARRAY ( 12r 50»75 ) » 

.  NNSAV(75)rSIG(12»75)»MYEAR(50) 

DATA  ARRAY/4500010./ 

C  *****  NEXT  SECTION  ALLOWS  FOR  MINIMUM  MONTHLY  STORAGES  TO  BE  PRINTED  OR 
C  *****  SKIPPED*  DEPENDING  ON  THE  CONTROL  VARIABLE  JMIN 
IF ( JMIN-1 )  510*500 
500  NFK=NOUT 
GO  TO  520 

510  NFK=NOUT-NRES 
520  CONTINUE 

C  *****  WRITE  THE  LAST  TWO  DAYS  FLOW  VALUES  ON  THE  BINARY  TAPE 
WRITE( NTAPE ) I YR2  * IM02* IDY2*(Q2(N) *N-1 *NOG ) 

WRITE ( NTAPE) I YR1 ,IM01*IDY1*(Q1(N)*N=1 *NOG ) 

IF(NYEAR.EQ .51 >  NYEAR-50 

C  *****  FIND  THE  SLOT  NUMBER  OF  THE  SFLECTED  OUTPUTS  ON  THE  BINARY  TAPE 
DO  20  N=1 *NOUT 
DO  10  NSAVE=1 *NOG 

IF(NOUTPT(N) . EQ.NAhEG(NSAVE) )  GO  TO  15 
10  CONTINUE 

WRITE! 6 • 12 ) 

12  FORMAT ( 44H  ERROR  IN  OUTPUT  S  SELECTED  OUTPUT  NOT  FOUND  ) 

15  CONTINUE 

NNSAV(N)=NSAVE 
20  CONTINUE 

REWIND  NTAPE 
READ(NTAPE)  NOG 

READ(NTAPE)  (NAMEG( I ) » Iel *NOG) 

21  CONTINUE 

READ(NTAPE)  I YR* IMO* IDY* (Q( II ) *  11*1 *NOG> 

IOLD=IMO 

IBEG=1901 

IF(IMO.EQ.l)  IBEG=1900 

C  *****  BEGIN  BUILDING  THE  MONTHLY  OUTPUT  ARRAY 
DO  60  K=1*NYFAR 
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MYEAR(K)=IBEG+IYR 
DO  60  NM0=1 *  12 
DO  22  J=1 *  NDU  T 
NX=NNSAV<  J) 

ARRAY(NMO»K. J)=G(NX ) 

22  CONTINUE 
DO  30  I"2»32 

READ (NT APE )  I YR. IHO. IDY , ( 0( 1 1 ) , II  =  1  .NOG) 

IF (EOF (NT APE )) 70.25 
25  CONTINUE 

C  ttttt  CHECK  TO  SEE  IF  A  NEW  MONTHS  RECORD  HAS  BEGUN 
IF(INO.NE.IOLD)  GO  TO  40 

C  *****  FILL  THE  OUTPUT  ARRAY  WITH  THE  SUM  OF  DAILY  VALUES 
DO  30  IN=1 .NOUT 
NX-NNSAV(IN) 

ARRAY(NMO»K» IN) = ARRAY (NMO.K. IN)+Q(NX> 

30  CONTINUE 

40  CONTINUE 

DAY=I-1 

C  *****  DIVIDE  THE  SUM  OF  THE  DAILY  FLOWS  BY  THE  NO.  OF  DAYS  IN  THE  MONTH 
DO  50  IX=1 »NOUT 

ARRAY ( NMO , K . I X ) = ARRAY (NMO.K. IX) /DAY 
C  *****  KEEP  TRACK  OF  MONTHLY  STATISTICS 

FHAX(NMO» IX )=AMAX1 ( ARRAY(NMO»K » IX ) »FMAX(NMO. IX) ) 

FNIN(NNO» IX)=AMIN1 ( ARRAY ( NMO.K. IX) .FMIN(NMO.IX)) 

SIG(NMO> IX )=SIG( NMO. IXH  ARRAY ( NMO.K. IX  >  KARR AY (NMO.K . IX ) 

FMEAN(NMO> IX)=FMEAN(NMO.IX)+ARRAY(NMO»K» IX) 

50  CONTINUE 

IOLD=IMO 
60  CONTINUE 

C  *****  FINISH  CALCULATIONS  FOR  LAST  DAY  ON  THE  BINARY  TAPE 
70  CONTINUE 

YEAR=NYEAR 

DAY-I-1 

DO  80  IY=1.N0UT 

ARRAY(NHO*K* I Y ) =ARRAY(NMQ» K> I Y ) /DAY 
FMAX(NM0.IY)=AMAX1< ARRAY (NMO.K. IY) .FMAXINMO. IY)) 

FMIN(NMO» I Y )=AMIN1 ( ARRAY (NMO.K? I Y) »FMIN(NMO» IY) ) 

SIG(NMO? I Y) =SIG(NMO» IY)tARRAY(NM0*K» I Y ) tARRAY ( NMO. K » I Y ) 

FMEAN(NMO. I Y )=FMEAN(NMO» IY ) + ARRAY (NMO.K. IY ) 

80  CONTINUE 

C  *****  DO  THE  STATISTICAL  CALCULATIONS 
DO  90  J=1 .NOUT 
DO  90  JJ=1 .12 

C  *****  CALCULATE  THE  MONTHLY  MEANS  OF  NYEAR  MONTHS  OF  FLOW 
FNEAN<JJ»J)=FMEAN<JJ.J)/YEAR 
FM*FMEAN< JJ. J)*FMEAN< JJ. J> 

C  *****  CALCULATE  THE  STANDARD  DEVIATIONS  OF  THE  MONTHLY  VALUES 
FS“SIG( JJ. J)/YEAR 
IF(FS-FM)  88.88.85 
SIG(JJ.J)»SQRT(FS-FM) 
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60  TO  90 

88  SIG<JJ»J>=0. 

90  CONTINUE 

C  *****  BEGIN  OUTPUT  OF  TABLES 
DO  140  J=1 »NOUT 

C  *****  WRITE  HEADING  OF  OUTPUT  TABLE 

WRITE <6.100 )(NTITLE3(JJ»J)»JJ= 1.10) 

100  FORMAT ( 1H1 .3X.30HCDM  /  WATER  RESOURCES  DIVISION. 19X. 

.  29H*  t***FINAL***t*  >20X. 

.  31HDELAWARE  RIVER  DAILY  FLOW  MODEL. /»4X» 

.  24HANNANDALE.  VIRGINIA  .64X.10A4) 

IF ( IBEG.EQ. 1901 )  GO  TO  109 
WRITEU.108) 

108  FORMAT (/.55X.23H  AVERAGE  MONTHLY  VALUES. /.4X.8HCALENDAR./. 

.  5X.4HYEAR.6X.3HJAN.7X.3HFEB.7X.3HMAR.7X.3HAPR.7X.3HMAY. 

.  6X.4HJUNE.6X.4HJULY.7X.3HAUG.7X.3HSEP.7X.3H0CT .7X.3HN0V. 

.  7X.3HDEC./) 

GO  TO  111 

109  CONTINUE 
WRITE<6.110) 

110  FORMAT (/.55X.23H  AVERAGE  MONTHLY  VALUES. /.4X.5HWATER./. 

.  5X . 4HYEAR . 6X . 3H0CT . 7X . 3HN0V . 7X . 3HDEC . 7X » 3H JAN . 7X . 3HFEB . 

.  7X , 3HMAR . 7X . 3HAPR . 7X . 3HMAY , 6X . 4H JUNE . AX , 4H JULY . 7X . 3HAUG . 

.  7X.3HSEP./) 

Ill  CONTINUE 

C  *****  WRITE  THE  MONTHLY  FLOW  ARRAY 
DO  130  JX=1 .NYEAR 
IF( J-NFK)320»320»330 
330  IFK=J-NFK 

DO  340  NM0=1 . 12 

340  ARRAY ( NMO . JX . J ) =QMM (NMO. JX.IFK) 

320  CONTINUE 

WRITE<6» 120)MYEAR( JX) * (ARRAY (NMO. JX.J) »NM0=1 .12) 

120  FORMAT <3X» 16. 12F10. 1 ) 

125  FORMATS, 6X.3HAVG.12F10.1) 

130  CONTINUE 

C  *****  WRITE  MONTHLY  AVERAGE  OF  NYEAR  VALUES 
WR I TE ( 6 . 1 25  > ( FMEAN ( NMO . J  > . NMO= 1 . 1 2 ) 

140  CONTINUE 

WRITER, 150) 

150  FORMAT (////, 27H  SIMULATION  ENDED  NORMALLY  ) 

RETURN 

END 
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Subroutine  OUTPUT 


Identical  to  Subroutine  MONTH  in  the  means  of  connecting  nodes  to  out¬ 
put  and  error  messages.  Subroutine  OUTPUT  is  different  in  that  it  will  print  out 
the  daily  values  for  each  requested  output.  Figure  11-11  shows  the  minor  dif¬ 
ferences  in  the  two  flowcharts.  The  tables  will  appear  beginning  with  the 
earliest  year  to  the  final  year  of  simulation,  completing  an  entire  selected 
output's  results  before  beginning  the  next.  Each  table  is  arranged  in  rows  by 
the  day  of  the  month  and  in  columns  by  month.  Also  calcualted  for  each  month 
are  the  mean  and  standard  deviation,  minimum  and  maximum  values,  and  the  monthly 
total . 

Subroutine  OUTPUT  ends  the  simulation  in  the  same  manner  as  MONTH:  . 
SIMULATION  ENDED  NORMALLY. 
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Figure  1 1  - 11 .  Flowchart  for  Subroutine  OUTPUT. 
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SUBROUTINE  OUTPUT 

ttmtmtmtmmmmmsmsmmsmmtmmtmtmmmm 

THE  OUTPUT  ROUTINE  PRODUCES  YEARLY  TABULATIONS  OF  DAILY  FLOWS* 
OUTPUT  REQUIRES  A  REWIND  OF  THE  UNFORMATTED  OUTPUT  TAPE* 

luutmmmmttmmmtummutmmmtmmmnmutm 

COMMON  /CONTROL/  NYEAR, NT  APE » NOUT » NOG  »NN0DE , ITAPE » IYR, IMO, IDY 
COMMON  /CHANL/  NGT0( 100) . NAMEG< 100) ,NGT0G(100, 10) ,N0UTPT( 100) , 

.  LINK(IOO) »Q(100) , GIN( 100) ,QFK(4) ,QFK0(4>  »QPR(4> » 

.  GPRG<4) 

COMMON  /INOUT/  NTITLEI(10),NTITLE2(10>,NTITLE3(10»100) 

DIMENSION  SIG( 12) »FMAX (12),FMIN(12) »FDAY (12) »FMEAN( 12) » 

.  FTOT< 12) » ARRAY (31 » 12) , OUT (8*100) ,NNSAV( 100) 

COMMON  /RESERU/  NRES,XFL0W,IFLAG,V0L(10> »LAGRES(10) ,LNKRES(10) . 

.  EVAP<12»10) ,ST0RE(10)»C0NS(12»10),Q01(100), 

.  002(100). NOGG . ICK . FLOHL . ZFLOU » ICK 1 

COMMON  /DATE/  I YR3.IM03.IDY3, IYR2.IM02. IDY2. IYR1 . IM01 » IDY1 . 

.  01(100.02(100) 

COMMON  /CAL/  CAL IB1 . JOUT, JMIN 

DIMENSION  KMO( 12) 

DATA  KMO  /31. 28. 31. 30. 31. 30. 31. 31. 30. 31. 30. 31/ 

DO  9000  N0=1 .NOUT 

C  *****  FIND  THE  LOCATION  OF  THE  SELECTED  OUTPUT  ON  THE  UNFORMATTED  TAPE 
DO  30  NSAVE=1 .NOG 

IF  (NOUTPT (NO) .EQ.NAMEG(NSAUF) )  GO  TO  35 
30  CONTINUE 
WRITE(4,33> 

33  FORMAT (17H  ERROR  IN  OUTPUT  ) 

35  CONTINUE 

NNSAV(NO)=NSAVE 
REWIND  NT APE 
READ (NT APE) NOG 

RE AD (NT APE) (NAMEG( I ) . 1=1 .NOG) 

READ(NTAPE) IYR. IMO. IDY . (Q( II ) > 11=1 .NOG) 

37  CONTINUE 

C  *****  LOOP  ON  NUMBER  OF  YEARS  OF  THE  RUN 
DO  8000  K=1 .NYEAR 
ARRAY(29.2)=0. 

DO  50  J=1 » 12 

IF( J.NE. IMO)GO  TO  45 

ARR AY ( 1 . J ) =Q ( NSAVE ) 

DO  40  1=2.32 

RE AD (NT APE) IYR. IMO, IDY .(Q(1I).II=1 ,NOG) 

IF  (EOF (NTAPE) )  38,39 

38  IM0=13 
IYR=IYRF1 

39  CONTINUE 

IF  (IMO.NE.J)  GO  TO  50 
ARRAY(I,J)*Q( NSAVE) 

40  CONTINUE 
GO  TO  50 

C  *****  ZERO  SUMMERS 
45  DO  47  1*1,31 
ARRAY(I» J)=0, 

47  CONTINUE 
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50  CONTINUE 
55  CONTINUE 

DO  100  J=l,12 
FMIN(J)*1.E30 
FNEAN(J>=0. 

FT0T< J)=0. 

FMAX(J>=0, 

FDAY(J)=0, 

SIG< J)=0. 

100  CONTINUE 

C  I**#*  PERFORM  STATISTICAL  CALCULATIONS 
IIYR*IYR-1 
DO  200  J=1 .12 
LL=KN0( J) 

IF( J.NE.2)G0  TO  110 
IF<MOD<IIYR,4>.EQ.O>LL=LLF1 
110  CONTINUE 

DO  150  Ial »LL 

120  FMEAN(J>aFMEAN< J)+ARRAY< I » J) 

SIG<J)=SIG(J) +ARRA Y ( I » J ) < ARRAY ( I  *  J ) 

FMAX ( J ) =AMAX1 <  ARRAY (ItJ) *FMAX( J) ) 

FNIN( J)aAMINl (ARRAY(I» J) ,FMIN( J) ) 

FDAY ( J)=FDAY ( J) +1 . 

150  CONTINUE 
200  CONTINUE 

NNAME=NOUTPT(NO) 

WRITE(6,7280)NTITLE1 ,NTITLE2,NNAME»  <NTITLE3(N»N0) »N=1 , 10) 

7280  FORMAT ( 1H1 » /»3X»30HCDM  /  HATER  RESOURCES  DIVISION, 20X, 

*  29Hf  mtFINAltm  *>9X»10A4»/»3X» 

*  21HANNANDALE i  VIRGINIA  ,29X, 

*  29H*  DELAWARE  RIVER  t,9X,10A4,/3X, 

*  7HN0DE  ,I5,38X,29H*  DAILY  FLOH  MODEL  *  ,9X, 

*  10A4) 

MM=l900+l YR-1 
URITE<6,7260)  MM 

7260  FORMAT ( /// » 24 ( 2H — ) » 35HAVERAGE  DAILY  VALUES  SIMULATED  FOR  ,14* 

.  23(2H— ),/) 

URITE(6,7310) 

7310  FORMAT ( / / , 3X, 3HDAY ,  FX , 3H JAN i 7X , 3HFED , 7X , 3HMAR , 7X , 3HAPR , 7X , 3HMAY , 

*  7X , 3H JUN , 7X , 3N JUL , 7 X , 3HAUG , 7X , 3HSEP , 7X , 3H0CT , 7X , 

t  3HN0V, 7X»3HDEC) 

WRITE(<5,73O0)  (NDAY,  (ARRAY  (MDAY ,MON) »M0N=1 , 12) ,MDAY*1 ,31 ) 

7300  FORMAT ( I5,3X,12F10t2) 

CALCULATE  DEVIATIONS 

DO  400  J»l,12 

IF  (FDAY(J).LE.O.O)  GO  TO  400 
FT=FMEAN< J) /FDAY( J) 

FTOT  ( J)aFMEAN( J) 

FMEAN( J)"FT 
FT*FTtFT 

FT1=SIG(J)/FDAY(J> 

IF  (FT1-FT)  378,378,379 
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379  SIG( J)=SQRT (FT1-FT) 

GO  TO  400 
378  SIG(J)=0. 

400  CONTINUE 

DO  450  1=1,12 

IF  <FHIN<I).GE.1.0E29>  FMIN(I)=0. 

450  CONTINUE 

WRITE <6, 7221) (FMIN(I), 1=1,12), (FMAX( J) » J=1 , 12) t 
.  (FMEAN(KO) »K0=1 ,12) 

7221  FORMAT </,5X,3HHIN,  12F10.2/,5X,3HMAX,  12F10.2, 

.  /4X»  4HMEAN, 12F 10 . 2 ) 

URITE(6,7222)  (FTOT(N) , N=1 , 12) , (SIG(L> ,L=1 , 12) 

7222  FORMAT ( 3X, 5HT0T AL » 12F 10.2, /3X,5HS  DEV,12F10.2) 
900  CONTINUE 

8000  CONTINUE 
9000  CONTINUE 

WRITE<6»9100) 

9100  FORMAT (26H1SIMULATI0N  ENDED  NORMALLY) 

STOP 

END 
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III.  INPUT  DESCRIPTION 


This  chapter  provides  a  detailed  description  of  the  card  data  require¬ 
ments  for  the  Daily  Model  along  with  an  example.  Table  III-l  gives  the  indivi¬ 
dual  card  structure,  including  variable  names  and  required  format.  This  program 
does  check  for  inclusion  and  order  of  the  card  groups.  Therefore,  an  iden¬ 
tifier,  A  through  M,  has  been  incorporated  into  each  data  block  in  the  first 
column. 


The  input  data  to  the  program  includes  the  card  data  deck  and  the 
Incremental  Inflow  Tape  ( I  IT) .  The  1 1 T  contains  the  average  daily  incremental 
inflow  to  nodes  in  the  entire  Delaware  River  Basin.  Table  1 1 1 -2  gives  the  loca¬ 
tion  and  model  node  numbers  along  with  the  corresponding  USGS  station  number. 
Some  of  the  model  nodes  do  not  have  USGS  stations  and  are  represented  by 
asterisks.  These  exact  model  node  numbers  must  be  used  in  the  card  input  deck 
or  LINKER  will  not  be  able  to  link  the  incremental  inflows  on  tape  to  the  spe¬ 
cified  model  node.  The  I  IT  is  the  driving  force  of  the  model  while  the  card 
input  deck  contains  program  control  information. 

Title  Cards  (A) 

The  first  two  cards  in  the  data  set  are  the  title  cards.  They  are 
printed  in  the  heading  of  the  daily  output  but  not  the  monthly.  They  contain  40 
alphanumeric  characters,  and  can  be  used  to  describe  the  particulars  of  the 
simulation  run. 
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TABLE  II I -1 
INPUT  DATA 
CARD  FORMAT 


Card  Card 

Group  No. 

Format 

Variable 

Name 

Description 

TITLE  CARDS 

A  1 

A1.4X 

AC 

Card  Identifier  (A) 

10A4 

NTITLE1 

Title  Card 

2 

A1.4X 

AC 

Card  Identifier  (A) 

10A4 

NTITLE2 

Title  Card 

CONTROL  CARDS 

B  1 

A1.4X 

BC 

Card  Identifier  (B) 

15 

I  TAPE 

Input  Tape  Number 

15 

NTAPE 

Output  Tape  Number 

15 

NYEAR 

Number  of  Years  to  be  Simulated 

F5.2 

CALI B I 

Multiplication  factor  for  the 
releases  from  NYC  reservoirs 

15 

NYCRES 

Number  of  New  York  City 
Reservoirs 

15 

JOUT 

Specifies  output  option  (-1 
daily  only,  0  monthly  only,  +1 
both) 

15 

JMIN 

Specifies  option  of  minimum 
monthly  reservoir  storages: 

+1  no  minimum  storages,  0  mini¬ 
mum  storages  in  special  output 

15 

MWL 

Wallenpaupack  option:  +1  use 
special  releases,  0  use  natural 
rel eases 

CHANNEL  CARDS 

C  1 

A1,4X 

CC 

Card  Identifier  (C) 

no 

NAMEG(I) 

Upstream  Node  Number 

no 

NGTO(I) 

Downstream  Node  Number 

E10.4 

XM(  I ) 

Lag  function  coefficient  de¬ 
fined  in  hours 

F10.0 

B  ( I ) 

Maximum  route  time  to  next  node 
in  hours 

no 

ICKK 

Flag  for  Objective  Node.  Equal 
to  1  if  is  1st  objective. 

Equal  to  2  if  is  2nd  objective. 
Otherwise,  blank. 

Card  1  repeated  for  each  node  to  be  modeled. 

To  specify  end  to  C  Cards,  one  additional  card  with  NAMEG(I)  equal  to  zero  Is 
required. 
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TABLE  II I -1 
INPUT  DATA 
CARD  FORMAT 
(Cont'd) 


Card  Card 

Group  No. 

Format 

Variable 

Name 

Description 

OUTPUT  NODE  CARDS 

D  1 

A1,4X 

DC 

Card  Identifier  (D) 

15 

NOUTPT(I) 

Output  Node  Number 

10A4 

NTITLE3(N, I ) 

Title  to  be  Printed 

Card 

1  repeated  for  each  nodal 

output  desired. 

To  specify  end  to  D 

Cards,  one 

:  additional  card  with  NOUTPT(I)  equal  to  zero  i: 

required. 

RESERVOIR  CARDS 

E  1 

A2,3X 

EC 

Card  Identifier  (El) 

15 

IRES 

Reservoir  Storage  Node  Number 

E10.4 

VOL ( I ) 

Maximum  Volume  (ft3) 

15 

LRES 

Node  Number  to  Which  IRES 

F10.2 

AREA 

Releases 

Surface  Area  Used  to  Calculate 

5A4 

NAMR ( 5 ) 

Evaporation  from  in/mo.  to 
cfs-day  (acres) 

Name  of  Reservoir  -  20 

2 

5X.12F5.0 

CONS (N, I ) 

Characters 

Conservation  Release  (cfs) 

3 

5X , 12F5.0 

ACONS(N.I) 

Monthly  Augmented  Conservation 

4 

5X.12F5.0 

E VAP ( N ,  I ) 

Release  (cfs) 

Evaporation  Rate  (in/mo) 

Cards  1,  2,  3, 

4  repeated  for 

each  reservoir. 

To  specify  end  to  E  cards,  one  additional  card  with  IRES  equal  to  zero  is 
required. 

If  NYCRES  equals  zero,  skip  F,  G,  H,  and  I  groups. 

SYSTEM  CONDITION  CARDS 


1 

A1.4X 

FC 

Card  Identifier  (F) 

15 

NCARDS 

Number  of  Cards  to  Follow 

2 

5X.I5 

NMO 

Month 

15 

NDY 

Day 

F10.0 

FLOW ( I ) 

Storage  Level  (bg)  that 

seperates  Normal  and  Drought 
Warning  Level s 


Repeat  Card  2  NCARDS  times. 
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TABLE  III-l 
INPUT  DATA 
CARD  FORMAT 
(Cont'd) 


Card 

Card 

Variable 

Description 

Group 

No. 

Format 

Name 

G 

I 

A1.4X 

GC 

Card  Identifier  (G) 

15 

NCARDS 

Number  of  Cards  to  Follow 

2 

5X  ,15 

NMO 

Month 

15 

NDY 

Day 

F10.0 

FLOW(I) 

Storage  Level  (bg)  that 
seperates  Drought  Warning  and 
Drought  Levels. 

Repeat  Card  2  NCARDS  times. 

DIVERSION  CARDS 

H 

1 

A1,4X 

CH 

Card  Identifier  (H) 

3F10.0 

DIVERS ( 3 ) 

New  York  City  Diversions  for  3 
Conditions  Normal,  Drought 
Warning,  Drought  Respectively 

(mgd) 

MONTAGUE 

OBJECTIVES  CARD 

I 

1 

A1.4X 

Cl 

Card  Identifier  (I) 

3FI0.0 

OBJECT ( 3 ) 

Montague  Flow  Objectives  for  3 
Conditions  Normal,  Drought 

Warning,  Drought  Respectively 
(cfs) 

SPECIAL 

OUTPUT  CARDS 

J 

1 

A1.4X 

CJ 

Card  Identifier  (J) 

15 

NSPO 

Number  of  Special  Outputs 

2 

A1.4X 

CJ 

Card  Identifier  (J) 

15 

NOUTPT(I) 

Special  Output  Flag  Number 

10A4 

NTITLE3(N,I) 

Description 

Card 

2  repeated  for  each  special 

output. 

D  &  R  CANAL  DIVERSION 

CARD 

K 

I 

A1,4X 

CK 

Card  Identifier  (K ) 

3F10.0 

DFLOW(3) 

D  &  R  Canal  Diversion  (cfs)  for 

Normal,  Drought  Warning  and 
Drought  conditions, 
Respectively 


If  NYCRES  equals  zero,  skip  L  group. 
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TABLE  III-l 
INPUT  DATA 
CARD  FORMAT 
(Cont*  d) 


Card 

Card 

Variable 

Description 

Group 

No. 

Format 

Name 

TRENTON 

OBJECTIVES  CARD 

L 

1 

A1,4X 

CL 

Card  Identifier  (L) 

3F10.0 

TARG ( 3 ) 

Trenton  flow  objectives  (cfs) 

for  Normal,  Drought  Warning  and 
Drought  conditions. 

Respectively 

CONSUMPTIVE  USES 

CARDS 

M 

1 

A1,4X 

CM 

Card  Identifier  (M) 

no 

NCON 

Number  of  nodes  with 
consumptive  uses  (losses) 

2 

A5 

CM 

Card  Identifier  (M) 

no 

ICON ( I ) 

Node  Number  for  subtraction  of 
consumptive  uses 

I2F5.0 

CLOSS ( I ,J ) 

Consumptive  use  for  each 

ICON ( I )  node,  one  for  each 
month  (cfs) 

Card  2  repeated  for  each  node  with  consumptive  use  data 
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TABLE  1 1 1-2 
MODEL  NODE  NUMBER 
AND 

USGS  STATION  LOCATION 


Model  Node 

USGS  Number 

Location 

1600 

01417000 

E  Branch  Delaware  River  at  Downsville,  NY* 

(1000)+ 

1010 

01421000 

E  Branch  Delaware  River  at  Fishs  Eddy,  NY 

1610 

01425000 

W  Branch  Delaware  River  at  Stilesvil le.NY? 

(1020)+ 

1030 

01426500 

W  Branch  Delaware  River  at  Hale  Eddy,  NY 

1035 

01427405 

Delaware  River  near  Cal  li coon,  NY 

1040 

01428500 

Delaware  River  near  Barryville,  NY 

1050 

01429000 

W  Branch  Lackawaxen  River  at  Prompton,  PA 

1060 

01429500 

Dyberry  Creek  near  Honesdale,  PA 

1070 

01430000 

Lackawaxen  River  near  Honesdale,  PA 

1080 

0143500 

Lackawaxen  River  at  Hawley,  PA 

1090 

01432000 

Wal lenpaupack  Creek  at  Wilsonville,  PA 

1100 

01434000 

Delaware  River  at  Port  Jervis,  NY 

1160 

01436000 

Neversink  River  at  Neversink,  NY  3 

( 1 1 10) + 

1120 

01437000 

Neversink  River  at  Oakland  Valley,  NY 

1130 

01438500 

Delaware  River  at  Montague,  NY 

1140 

01440200 

Delaware  River  below  Tacks  Island  Damsite 

1150 

01446500 

Delaware  River  at  Belvidere,  NJ 

^Pepacton  Reservoir  +If  the  NYC  Reservoirs  are  not  to  be  used  as 

2Cannonsvil le  Reservoir  reservoirs,  then  the  node  number  in  paren- 

^Neversink  Reservoir  thesis  must  be  used. 
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TABLE  1 1 1-2 
MODEL  NODE  NUMBER 
AND 

USGS  STATION  LOCATION 
(Cont'd) 


Model  Node 

USGS  Number 

Location 

1160 

01447800 

Lehigh  River  near  White  Haven,  PA 

1170 

01449800 

Pohopoco  Creek  below  Beltzville  Damsite,  PA 

1180 

* 

Aquasicola  Creek  at  Aquashicola  Damsite,  PA 

1190 

01450500 

Aquashicola  Creek  at  Palmerton,  PA 

1200 

01451000 

Lehigh  River  at  Walnutport,  PA 

1210 

01451800 

Jordan  Creek  near  Schnecksvi 1 le,  PA 

1220 

01452000 

Jordan  Creek  at  Allentown,  PA 

1230 

01453000 

Lehigh  River  at  Bethlehem,  PA 

1240 

01454700 

Lehigh  River  at  Glendon,  PA 

1250 

01456000 

Musconetcong  River  near  Hackettstown,  NJ 

1260 

01457500 

Delaware  River  at  Riegel svil le,  NJ 

1280 

01459500 

Tohickon  Creek  near  Pipersville,  PA 

1290 

* 

Delaware  River  near  Point  Pleasant 

1300 

01460500 

D  &  R  Canal  at  Kingston,  NJ 

1310 

01463500 

Delaware  River  at  Trenton,  NJ 

1320 

01464000 

Assunpink  Creek  at  Trenton,  NJ 

1330 

01464500 

Crosswicks  Creek  at  Extonville,  NJ 

1340 

01465500 

Nashaminy  Creek  near  Langhorne,  PA 

1350 

01467000 

N  Branch  Rancocas  Creek  at  Pemberton,  NJ 

1360 

★ 

Delaware  River  near  Torresdale,  PA 

1370 

* 

Delaware  River  near  Pier  II  North 

♦Not  a  USGS  Station 
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TABLE  1 1 1-2 
MODEL  NODE  NUMBER 
AND 

USGS  STATION  LOCATION 
(Cont'd) 


Model  Node 

USGS  Number 

Location 

1380 

01467500 

Schuylkill  River  at  Pottsville,  PA 

1390 

01467950 

Schuylkill  River  at  Cressona.PA 

1400 

01468500 

Schuylkill  River  at  Landi ngvil le,  PA 

1410 

01469500 

Lower  Schuylkill  River  at  Tamaqua,  PA 

1420 

01470000 

Lower  Schuylkill  River  at  Drehersvil le,  PA 

1430 

01470500 

Schuylkill  River  at  Berne,  PA 

1440 

01470756 

Maiden  Creek  at  Virginville,  PA 

1450 

01470960 

Tulpehocken  Creek  at  Blue  Marsh,  PA 

1460 

01471000 

Tulpehocken  Creek  near  Reading,  PA 

1470 

01471500 

Schuylkill  River  near  Reading,  PA 

1480 

01472000 

Schuylkill  River  at  Puttstown,  PA 

1490 

★ 

Schuylkill  River  near  Limerick,  PA 

1500 

014723000 

Perk i omen  Creek  at  Graterford,  PA 

1510 

01474500 

Schuylkill  River  at  Philadelphia,  PA 

1520 

★ 

Delaware  River  below  Schuylkill  Confluence 

1530 

01477000 

Chester  Creek  near  Chester,  PA 

1540 

★ 

Delaware  River  near  Chester,  PA 

1550 

01479000 

White  Clay  Creek  near  Newark,  DE 

1560 

01480000 

Red  Clay  Creek  at  Wooddale,  DE 

1570 

01481000 

Brandywine  Creek  at  Chadds  Ford,  PA 

1580 

01481500 

Brandywine  Creek  at  Wilmington,  DE 

1590 

★ 

Delaware  River  at  Delaware  Memorial  Bridge 

*Not  a  USGS  Station 
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Control  Card  (B) 


This  third  card  contains  the  controlling  variables;  the  input  and  out¬ 
put  tape  numbers,  the  numbers  of  years  of  simulation,  the  multiplication  factor 
for  the  New  York  City  reservoirs,  the  number  of  New  York  City  reservoirs  and, 
the  option  of  monthly  or  daily  final  output.  A  negative  one  (-1)  indicates 
daily  output  only,  a  zero  (0)  specifies  monthly  output  only,  and  a  positive  (+1) 
gives  both  daily  and  monthly  output.  When  only  the  monthly  output  has  been  spe¬ 
cified,  another  output  option  is  available.  That  is  the  inclusion  of  minimum 
monthly  storages.  A  one  (1)  indicates  no  minimum  storages  are  requested  while  a 
zero  (0)  indicates  that  the  minimum  monthly  storages  for  every  reservoir  will  be 
included  in  the  Special  Outputs.  The  final  variable  controls  the  operation  of 
Wal lenpaupack.  Special  releases,  as  described  more  fully  in  the  Daily  Flow 
Model  Report^,  can  be  released  from  Wal lenpaupack  to  help  augment  the  Montague 
flow.  If  this  variable  is  a  one  (1),  the  special  releases  will  be  made.  A  zero 
(0)  indicates  the  natural  releases  from  Wal lenpaupack  are  desired. 

Channel  Cards  (C) 

These  are  the  node  input  cards,  specifying  the  channel  between  a  pair 
of  nodes,  one  upstream,  one  downstream.  One  card  is  required  for  each  node. 

The  user  should  note  that  each  upstream  node  is  connected  to  only  one  downstream 
node,  i.e.,  NAMEG(I)  is  not  duplicated.  However,  each  downstream  node  can  have 
multiple  upstream  nodes,  i.e.,  NGTO(I)  can  be  duplicated.  When  inputting  the 
node  numbers,  it  is  essential  that  they  be  put  in  order  from  upstream  to 
downstream.  The  first  flow  objective  node  connected  with  backup  is  flagged  in 
this  card  series  by  putting  a  one  (1)  in  column  35.  The  second  flow  objected 
node  connected  with  the  next  day  release  scheme  is  flagged  by  putting  a  two  (2) 
in  column  35.  A  maximum  number  of  100  model  nodes  are  possible.  Of  these  100 
nodes,  only  10  upstream  nodes  may  be  routed  into  any  one  downstream  node. 


^Development  of  A  Daily  Flow  Model  of  the  Delaware  River  Which 
Incorporates  Reservoir  Systems  Analysis,  Volume  II,  Camp  Dresser 
&  McKee,  March  19817 


1 1 1-9 


Output  Nodes  (D) 


Two  types  of  output  tables  can  be  generated.  One  is  the  actual  flows 
at  the  system's  nodes,  which  are  specified  by  the  D  cards.  The  other  is  a  spe¬ 
cial  output  described  by  the  J  cards.  The  D  card  contains  the  card  identifier, 
the  model  node  number  where  the  flow  is  requested  plus  an  alphanumeric  descrip¬ 
tion  of  the  node.  A  maximum  of  75  output  tables  can  be  generated.  This  inclu¬ 
des  both  the  Output  Node  tables  plus  the  Special  Output  tables  (D  cards  and  J 
cards) . 

Reservoir  Input  (E) 

This  group  contains  a  series  of  four  cards  for  each  reservoir.  If  no 
reservoirs  are  to  be  used,  two  items  must  be  changed.  Only  one  reservoir  card 
is  input:  an  "E"  followed  by  a  blank.  Also,  the  NYC  reservoirs’  node  numbers 
must  be  replaced  by  the  node  numbers  in  parentheses  in  Table  1 1 1-2.  When  reser¬ 
voirs  are  included,  the  first  card  contains  the  reservoir  storage  node  number, 
the  maximum  volume  in  cubic  feet,  the  node  to  where  the  release  is  made  (usually 
the  inflow  node),  and  the  surface  area  used  to  calculate  the  evaporation  in  cfs- 
day.  The  name  of  the  reservoir  is  included  to  help  with  bookkeeping.  It  is 
essential  to  specify  the  reservoir  storage  node  as  1000  plus  the  number  of  the 
node  where  the  incremental  inflow  is  associated.  This  inflow  node  is  generally 
equal  to  the  release  node.  A  maximum  of  ten  reservoirs  can  concurrently  be  used 
in  the  model.  The  second  card  contains  the  average  daily  basic  conservation 
releases  in  cfs  for  each  month.  The  third  card  contains  the  average  daily 
augmented  conservation  release  in  cfs.  The  fourth  card  contains  the  average 
daily  evaporation  rates  in  inches/month  for  each  month. 
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The  first  three  reservoirs  must  always  be  the  New  York  City  reservoirs, 
Pepacton,  Neversink  and  Cannonsville  unless  zero  NYC  reservoirs  have  been  spe¬ 
cified  on  the  B  Card.  If  more  than  three  reservoirs  are  used  in  the  simulation, 
the  New  York  City  reservoirs  must  be  placed  first. 

The  special  cases  of  Merrill  Creek  and  Blue  Marsh  are  also  coded  on  the 
E  Cards.  These  two  reservoirs  are  not  operated  to  maintain  any  particular 
objective.  They  do,  however,  act  as  reservoirs  and  are  therefore  specified  in 
this  card  series.  If  Merrill  Creek  is  in  operation,  it  must  be  placed  before 
all  other  reservoirs  used  to  maintain  the  Trenton  objective  and  after  any  NYC 
reservoirs.  If  Blue  Marsh  is  on  line,  it  must  be  placed  after  all  other  reser¬ 
voirs  have  been  specified. 

System  Conditions:  New  York  City  Reservoir  Storage  Levels  (F  and  G) 

The  F  and  G  card  series  set  the  system’s  conditions  by  the  New  York 
City  reservoir  storage  levels.  If  there  are  no  NYC  reservoirs,  these  two  groups 
are  skipped.  The  first  (F  series)  designates  the  storage  levels  at  which  the 
system  changes  from  Normal  to  Drought  Warning  conditions. 

The  second  (G  series)  designates  the  storage  levels  at  which  the  system 
changes  from  Drought  Warning  to  Drought  conditions.  At  the  beginning  of  each 
day  during  simulation,  the  total  New  York  City  storage  in  the  system  is  compared 
to  the  reservoir  levels  specified  in  the  F  and  G  cards  for  that  particular  day. 
The  system  condition  is  then  set. 

The  first  card  in  each  set  has  the  identifier  and  the  number  of  cards 
to  follow.  The  next  n-cards  give  the  storage  levels  at  appropriate  dates.  The 
month,  day  and  the  total  storage  level  in  billion  gallons  are  given.  With  a 
beginning  point  of  January  1,  the  storage  levels  thereafter  are  specified  at  the 
inflection  points  seen  in  Figure  1-1.  December  31  is  specified  to  end  the  year. 
The  levels  at  each  day  are  calculated  by  linear  interpolation  in  the  program. 
From  these  daily  values,  the  System  Conditions  are  set  every  day. 
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Diversions  (H) 


As  part  of  the  New  York  City  reservoir  operation  matrix,  the  diversions 
to  New  York  City  are  specified  on  this  card.  Three  diversions  for  the  three 
reservoir  conditions.  Normal,  Drought  Warning  and  Drought,  are  specified  respec¬ 
tively.  It  should  be  noted  that  only  Pepacton,  Cannonsville  and  Neversink  have 
diverisons  taken  from  them.  Any  other  reservoir  contributes  only  to  the  Trenton 
flow  objective.  The  diversions  are  the  total  taken  from  the  system,  not  from 
each  reservoir.  The  units  are  in  mgd.  If  there  are  no  NYC  reservoirs,  this 
card  is  eliminated. 

Montague  Objectives  (I) 

The  Montague  target  flow  objective  (or  the  objective  connected  to 
backup)  is  specified  on  the  I  card.  Three  objectives  are  specified  in  cfs, 
again  relating  to  the  New  York  City  System  Conditions.  This  flow  objective  is 
used  to  determine  the  deficit  of  the  Montague  flow  and  thus  the  release  needed 
to  maintain  the  Montague  objective.  A  card  identifier  followed  by  the  three 
objectives  are  on  this  card.  As  with  the  H  card,  if  there  are  no  NYC  reser¬ 
voirs,  this  card  is  eliminated  from  the  input  deck. 

Special  Outputs  (J) 

The  Special  Outputs  have  been  included  in  the  model  to  give  information 
additional  to  the  flows  generated  by  the  model.  Table  1 1 1-3  lists  the  Special 
Outputs  currently  available  along  with  their  flag  number.  One  card  is  required 
for  each  Special  Output.  The  flag  number  must  be  specified  higher  than  the 
reservoir  storage  numbers  so  as  not  to  confuse  them. 

When  setting  these  up  in  the  input  card  deck,  they  must  be  in  the  same 
order  as  in  Table  1 1 1-3.  This  is  due  to  the  order  of  definition  of  each  of  the 
Special  Outputs  in  RESOP.  If  any  change  is  desired,  the  FORTRAN  of  the  program 
would  have  to  be  altered. 
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The  output  from  the  D  cards  is  in  cfs.  The  Special  Outputs  are  in 
several  units:  bg,  cfs,  and  mgd.  The  user  should  take  care  when  analyzing  the 
results  to  remember  the  different  units. 

TABLE  1 1 1-3 
SPECIAL  OUTPUTS 


Flag 

Description 

6000 

Releases  to  Maintain  Montague  Flow  (cfs) 

6010 

Releases  to  Maintain  Trenton  Flow  (cfs) 

6020 

Total  Diversions  to  NYC  (mgd) 

6030 

Combined  NYC  Reservoir  Storage  (bg) 

6040 

Combined  Lower  Basin  Reservoir  Storage  (bg) 
(this  includes  Prompton  and  excludes  Merrill 
Creek  and  Blue  Marsh) 

D  &  R  Canal  Diversions  (K) 

The  Delaware  and  Raritan  (D  &  R)  Canal  is  now  set  to  be  operated  by 
the  same  conditions  as  the  New  York  City  reservoirs.  Therefore,  three  values 
are  input  corresponding  to  diversions  during  Normal,  Drought  Warning  and  Drought 
Conditions.  These  diversions  are  specified  in  cfs.  If  there  are  no  NYC  reser¬ 
voirs,  the  K  group  is  not  included. 

Trenton  Objectives  (L) 

The  Trenton,  or  second  target  flow  objective  is  specified  on  the  L 
card.  Three  objectives  are  specified,  relating  to  the  New  York  City  System 
Conditions.  However,  because  the  Trenton  objective  is  not  governed  by  the  New 
York  City  reservoirs,  they  are  specified  the  same.  The  objectives  are  specified 
in  cfs.  Currently  set  at  Trenton,  this  second  objective  node  can  be  changed  by 
setting  ICKK  in  the  C  cards  equal  to  two  (2)  for  a  different  node. 


Consumptive  Uses  (M) 


As  a  prediction  of  future  conditions,  consumptive  uses  have  been  added 
to  the  simulation.  The  card  identifier  and  the  number  of  cards  to  follow  is 
specified  on  the  first  card.  The  next  n-cards  contain  the  node  number  (taken 
from  Table  II 1-2)  which  is  specified  first,  followed  by  average  daily  consump¬ 
tive  uses,  from  January  to  December,  in  cfs.  Each  day,  these  consumptive  uses 
will  be  subtracted  from  the  appropriate  node. 

Example  Input 

The  example  input  deck  which  follows  shows  the  data  described  in  this 
chapter.  The  results  are  given  in  Chapter  IV.  The  first  two  cards  are  the  A 
Cards  or  title  cards.  The  third  card  is  the  B  Card  or  the  control  card.  On 

this  card  is  specified  the  number  of  years  of  the  simulation  as  well  as  other 

flags  and  variables.  In  this  example,  there  are  four  years  in  the  simulation 
and  all  three  NYC  reservoirs  are  used.  Monthly  output  only  for  the  tabular 
results  is  flagged.  The  next  card  series  is  the  C  Cards  which  specify  all  of 
the  nodes  to  be  used  in  the  network.  Each  card  contains  the  upstream  node,  the 
downstream  node  and  the  two  variables  used  in  the  routing  function.  In  this 
case,  all  62  nodes  are  being  used.  Also  identified  are  the  target  nodes. 
Montague  (node  number  1130)  is  the  first  target  node.  Trenton  (node  number 
1310)  is  the  second  target  node.  The  C  Cards  are  terminated  by  a  blank  C  Card. 

Only  three  nodes'  flows  (D  Cards)  are  requested  for  tabular  output: 

Montague,  Trenton,  and  Delaware  Memorial  Bridge.  The  D  Cards  are  terminated  by 
a  blank  D  Card.  After  the  D  Cards  the  reservoir  information  (E  Cards)  is  spe¬ 
cified.  Each  reservoir  has  four  cards.  The  first  contains  the  storage  node 
number  along  with  surface  area  and  maximum  volume.  The  second  card  gives  the 
basic  conservation  releases  for  each  month  of  the  year.  The  third  card  gives 
the  augmented  conservation  releases  for  each  month  of  the  year.  The  final  card 
gives  the  evaporation  of  each  reservoir  in  inches  per  month  for  each  month  of 
the  year.  After  the  reservoir  information,  six  points  are  used  to  delineate 
each  of  the  System  Conditions  (F  and  G  Cards).  The  first  two  values  give  the 
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month  and  day  of  the  occurrence  of  the  last  value,  the  total  NYC  reservoir 
storage.  The  next  two  cards  set  the  diversions  to  NYC  (H  Card)  and  the  Montague 
objectives  (I  Card).  Each  of  these  is  linked  to  the  System  Conditions  and 
therefore  all  three  values  are  different.  The  next  card  series  is  the  special 
outputs  (J  Card).  Six  have  been  specified  here  to  give  a  clear  illustration  of 
the  different  kinds  of  possible  special  outputs.  The  D&R  Canal  diversions  (K 
Card)  are  specified  next.  These,  too,  are  linked  to  the  System  Conditions  and 
are  therefore  different.  Trenton's  objectives  (L  Card)  on  the  other  hand  are 
identical,  being  independent  of  the  System  Conditions.  Finally,  consumptive 
uses  (M  Card)  occur  at  seven  nodes.  The  node  number  is  followed  by  the  consump¬ 
tive  use  for  each  month  of  the  year.  Currently,  these  values  are  constant  for 
each  node. 


A  DOCUMENTAION  SIMULATION  EXAMPLE 
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IV.  OUTPUT  DESCRIPTION 


The  output  for  each  simulation  is  divided  into  three  major  sections. 

The  first  is  an  echo  of  the  input  deck,  interspersed  with  calculations  done 
within  the  INDATA  subroutine.  The  second  reflects  the  simulation  as  it  pro¬ 
ceeds.  In  this  second  section,  a  message  is  printed  every  day  the  Montague  flow 
was  below  the  specified  target.  The  third  output  section  is  the  tabular  results 
specified  by  the  D  and  J  Cards.  This  is  printed  after  the  entire  simulation  is 
finished. 

INPUT  ECHO 

The  subroutine  INDATA  echoes  all  input  data.  It  also  performs  a  few 
operations  by  data  manipulation  or  by  calling  another  subroutine.  In  the 
example  output  given  at  the  end  of  this  chapter,  the  A,  B,  and  C  Cards  are 
quickly  identified.  Immediately  following  the  C  Cards  are  the  results  from  the 
subroutine  ORDER. 

The  purpose  of  this  matrix  is  to  present  the  slot  numbers  which  flow 
into  every  node.  The  first  column  is  the  slot  number  associated  with  the  node 
specified  in  the  second  column.  The  third  column  is  the  node  into  which  the 
second  column  flows.  The  columns  afterwards  specify  the  slot  numbers  flowing 
into  the  secor/  column  node.  For  example,  the  12th  slot,  node  number  1100,  has 
three  nodes  flowing  into  it,  slot  numbers  6,  10  and  11  (node  numbers  1040,  1080 
and  1090  respectively).  After  the  matrix  is  printed,  ORDER  returns  to  INDATA, 
where  the  reading  of  the  input  data  is  resumed. 

The  D  Cards  are  printed  next.  In  this  simple  example  there  are  only 
three  nodes  whose  flows  will  be  printed:  Montague,  Trenton  and  the  Delaware 
Memorial  Bridge. 
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The  E  Cards,  reservoir  input,  show  some  minor  data  manipulation.  The 
evaporation  for  each  reservoir  is  changed  from  the  input  of  inches/month  to  cfs- 
day  using  the  surface  area  in  acres.  The  final  calculation  in  cfs-day  is 
printed  as  the  fifth  cards  in  the  E  Card  section  of  each  reseroivr. 

The  next  card  sections  F,  G,  H  and  I  Cards,  are  simple  input  echoes. 

The  J  Cards,  special  outputs,  are  those  described  in  Chapter  III.  Six  are  shown 
for  this  example  simulation.  The  K  and  L  Cards  are  simply  input  echoes  also. 
Seven  consumptive  uses  (M  Cards)  are  specified  in  this  example.  Generated  here 
are  the  slot  numbers  associated  with  the  consumptive  uses  rather  than  the  node 
number  which  was  input. 

The  next  three  lines  are  warnings,  specifying  those  nodes  on  the 
Incremental  Inflow  Tape  which  were  not  connected  to  any  model  node.  These  three 
"Warning"  nodes  are  the  naturalized  flows  for  the  NYC  reservoirs.  They  do  not 
include  the  special  net  inflow  data  for  the  reservoirs  which  are  necessary  when 
the  reservoirs  are  in  operation.  Therefore,  these  three  "Warning"  nodes  are 
specified  whenever  the  NYC  reservoirs  are  not  in  operation. 

The  next  seven  lines  are  the  dividing  line  between  the  first  section  of 
data  echo  and  the  second  section  of  actual  simulation.  These  seven  lines  are 
generated  in  LINKER  and  merely  reflect  whether  or  not  the  special  nodes  have 
been  specified.  These  special  nodes  are  described  fully  in  Chapter  III. 

ON-GOING  SIMULATION  OUTPUT 

This  next  section  of  the  example  shows  only  a  small  portion  of  the 
actual  output.  Each  time  Montague's  flow  is  less  than  the  target  and  the  model 
is  required  to  backup,  this  flag  is  printed  -  LOW  FLOW  AT  MONTAGUE  -  followed  by 
the  year,  month  and  day.  It  is  rather  monotonous  and  time  consuming  and  can  be 
easily  "switched  off"  within  the  program  by  simply  deleting  or  commenting  out 
the  write  statement. 


IV-2 


TABULAR  RESULTS 


Two  types  of  tables  can  be  generated  by  the  model.  A  monthly  table 
which  uses  the  daily  results  to  calculate  a  monthly  average  is  requested  in  the 
example  problem  presented  in  this  report.  A  daily  table  can  also  be  generated. 
Each  table  is  generated  for  each  output  requested  on  the  D  and  J  Cards. 

As  described  in  the  example  input,  only  monthly  average  results  for 
eight  tables  have  been  specified,  three  nodes'  flows  and  five  Special  Outputs. 
Two  have  been  put  on  a  single  page  for  ease  of  comparison.  In  reality,  each  is 
printed  on  a  separate  computer  page. 

Each  monthly  table  title  contains  a  brief  description  of  the  data  being 
presented.  The  tables  for  flows  will  have  the  USGS  station  number  and  name 
printed.  These  are  specified  in  the  D  Cards.  The  special  outputs  as  specified 
by  the  J  Cards  will  change  as  necessary.  It  is  recommended  to  always  specify 
the  units. 

The  table  is  arranged  by  water  year,  from  October  to  September.  Each 
month  in  each  year  has  the  monthly  average  value  printed.  At  the  end  of  the 
table,  the  average  value  over  all  the  years  for  each  month  is  calculated. 

An  example  of  the  daily  output  table  is  shown  following  the  monthly 
example  for  a  similar  simulation.  The  tables  are  slightly  different.  The  daily 
table's  title  contains  the  titles  from  the  A  Cards  as  well  as  those  given  on  the 
output  request  cards,  D  and  J.  This  table  contains  the  results  from  a  single 
year,  this  time  organized  by  calendar  year.  For  the  daily  output,  the  minimum, 
maximum,  total  and  standard  deviation  are  calculated  along  with  the  mean  value. 
These  are  not  included  in  the  monthly  tables  in  order  to  keep  a  50-year 
simulation's  results  on  a  single  page. 
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